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Abstract: 

In this work, we discuss the resonance states of a quantum particle in a periodic potential plus 
a static force. Originally this problem was formulated for a crystal electron subject to a static 
electric field and it is nowadays known as the Wannier-Stark problem. We describe a novel 
approach to the Wannier-Stark problem developed in recent years. This approach allows to 
compute the complex energy spectrum of a Wannier-Stark system as the poles of a rigorously 
constructed scattering matrix and solves the Wannier-Stark problem without any approximation. 
The suggested method is very efficient from the numerical point of view and has proven to be 
a powerful analytic tool for Wannier-Stark resonances appearing in different physical systems 
such as optical lattices or semiconductor superlattices. 
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Chapter 1 

Introduction 



The problem of a Bloch particle in the presence of additional external fields is as old as the 
quantum theory of solids. Nevertheless, the topics introduced in the early studies of the system, 
Bloch oscillations [||, Zener tunneling Q and the Wannier-Stark ladder are still the subject 
of current research. The literature on the field is vast and manifold, with different, sometimes 
unconnected lines of evolution. In this introduction we try to give a survey of the field, summarize 
the different theoretical approaches and discuss the experimental realizations of the system. It 
should be noted from the very beginning that most of the literature deals with one-dimensional 
single-particle descriptions of the system, which, however, capture the essential physics of real 
systems. Indeed, we will also work in this context. 



1.1. Wannier-Stark problem 

In the one-dimensional case the Hamiltonian of a Bloch particle in an additional external field, 
in the following referred to as the Wannier-Stark Hamiltonian, has the form 
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Hw = ^ + Vix)+Fx, V{x + d) = V{x), (1.1) 

where F stands for the static force induced by the external field. Clearly, the external field 
destroys the translational symmetry of the field-free Hamiltonian Hq = /2m + V{x). Instead, 
from an arbitrary eigenstate with Hy/^ = Eq^ , one can by a translation over / periods d 
construct a whole ladder of eigenstates with energies Ei = Eq + IdF, the so-called Wannier-Stark 
ladder. Any superposition of these states has an oscillatory evolution with the time period 

Tb = ^ . (1.2) 

known as the Bloch period. There has been a long-standing controversy about the existence 
of the Wannier-Stark ladder and Bloch oscillations P~p!g|], and only recently agreement about 
the nature of the Wannier-Stark ladder was reached. The history of this discussion is carefully 
summarized in [12,2^22|. 



Prom today's point of view the discussion mainly dealt with the effect of the single band 
approximation (effectively a projection on a subspace of the Hilbert space) on the spectral 
properties of the Wannier-Stark Hamiltonian. Within the single band approximation, the a'th 
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Figure 1.1. Schematic illustration of the Wannier-Stark ladder of resonances. The width of the 
levels is symbolized by the different strength of the lines. 



band of the field-free Hamiltonian Hq forms, if the field is applied, the Wannier-Stark ladder 
with the quantized energies 

E^^i = ea + dFl, l = 0±l,... , (1.3) 



where is the mean energy of the a-th band (see Sec. 1.2). This Wannier-Stark quantization 
was the main point to be disputed in the discussions mentioned above. The process, which is 
neglected in the single band approximation and which couples the bands, is Zener tunneling [Q]. 
For smooth potentials V{x), the band gap decreases with increasing band index. Hence, as the 
tunneling rate increases with decreasing band gap, the Bloch particles asymmetrically tend to 
tunnel to higher bands and the band population depletes with time (see Sec. p!^ ). This already 
gives a hint that Eq. (|1.3D can be only an approximation to the actual spectrum of the sytem. 



Indeed, it has been proven that the spectrum of the Hamiltonian ( |l.lD is continuous |23,24| 



Thus the discrete spectrum (1.3) can refer only to resonances p5|-|29|, and Eq. (1.3) should be 
corrected as 

£a,i = Eo^ + dFl-i^ , (1.4) 



(see Fig. LI). The eigenstates of the Hamiltonian (|l.lD corresponding to these complex energies, 



referred in what follows as the Wannier-Stark states ^a^i{x), are metastable states with the 
lifetime given by r = h/T^. To find the complex spectrum ( |1.4| ) (and corresponding eigenstates) 
is an ultimate aim of the Wannier-Stark problem. 
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Several attempts have been made to calculate the Wannier-Stark ladder of resonances. Some 
analytical results have been obtained for nonlocal potentials and for potentials with 



a finite number of gaps ||32H38|. (We note, however, that almost all periodic potentials have 
an infinite number of gaps.) A common numerical approach is the formalism of a transfer 
matrix to potentials which consist of piecewise constant or linear parts, eventually separated by 
delta function barriers |3^43]. Other methods approximate the periodic system by a finite one 



44-47]. Most of the results concerning Wannier-Stark systems, however, have been deduced 
from single- or finite-band approximations and strongly related tight-binding models. The main 
advantage of these models is that they, as well in the case of static (dc) field as in the cases 
of oscillatory (ac) and dc-ac fields [^9|-^, allow analytical solutions. Tight-binding models have 
been additionally used to investigate the effect of disorder [^7|-52|, noise [S3| or alternating site 



energies [|64|-|68| on the dynamics of Bloch particles in external fields. In two-band descriptions 



Zener tunneling has been studied |69-73], which leads to Rabi oscillations between Bloch bands 



1 74]. Because of the importance of tight-binding and single-band models for understanding the 



properties of Wannier-Stark resonances we shall discuss them in some more detail. 

1.2. Tight-binding model 

In a simple way, the tight-binding model can be introduced by using the so-called Wannier 
states (not to be confused with Wannier-Stark states), which are defined as follows. In the 
absence of a static field, the eigenstates of the field-free Hamiltonian, 

Ho = ^ + V{x)^ (1.5) 
are known to be the Bloch waves 

(t^aA^) = exp(iKx)Xa,«;(2;) , Xa^^ + d) = XaA^) ' (1-6) 

with the quasimomentum k defined in the first Brillouin zone — Tr/d < n < n/d. The functions 
(|1.6|) solve the eigenvalue equation 



Ho(j)a,Kix) = eaiK)(pa,^zix) , e„(K + 2-n/d) = eain) , (1.7) 

where €„(/«) are the Bloch bands. Without affecting the energy spectrum, the free phase of the 
Bloch function (pa, nix) can be chosen such that it is an analytic and periodic function of the 
quasimomentum k ]|7^. Then we can expand it in a Fourier series in k, where the expansion 
coefficients 

/Tv/d 
dKexp{-iKld)(t)a,K{x) (1.8) 
-n/d 

are the Wannier functions. 

Let us briefly recall the main properties of the Wannier and Bloch states. Both form orthogonal 
sets with respect to both indices. The Bloch functions are, in general, complex while the Wannier 
functions can be chosen to be real. While the Bloch states are extended over the whole coordinate 
space, the Wannier states are exponentially localized ]7G,77], essentially within the l-th. cell of 
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Figure 1.2. Left panel - lowest energy bands ea(K) for the potential V{x) = cos(3;) with parame- 
ters h = 1 and m = 1. Right panel - associated Wannier states ^o,o (solid line) and ipifi (dotted 
line). 



the potential. Furthermore, the Bloch functions are the eigenstates of the translation (over a 
lattice period) operator while the Wannier states satisfy the relation 

tpa,l+l{x) = 1pa,l{x - d) , (1.9) 



which directly follows from Eq. ( |1.8D . Finally, the Bloch states are eigenstates of Hq but the 
Wannier states are not. As an example. Fig. |1.2| shows the Bloch band spectrum ea{i^) and 
two Wannier functions ilJa,o{x) of the system ( |1.5D with V{x) = cosx, m = 1 and h = 1. The 
exponential decrease of the ground state is very fast, i.e. the relative occupancy of the adjacent 
wells is less than 10~^. For the second excited Wannier state it is a few percent. 

The localization property of the Wannier states suggests to use them as a basis for calcu- 
lating the matrix elements of the Wannier-Stark Hamiltonian ( [1.11 ). (Note that the field- free 
Hamiltonian ( |1.5| ) is diagonal in the band index a.) The tight-binding Hamiltonian is deduced 
in the following way. Considering a particular band a, one takes into account only the main 
and the first diagonals of the Hamiltonian Hq. From the field term x only the diagonal part is 
taken into account. Then, denoting the Wannier states resulting from the a-th band by |/), the 
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tight-binding Hamiltonian reads 

Htb = Y^i^a + dFl) \l){l\ + ^ ( |/ + + |/)(/ + 1| ) . (1.10) 



The Hamiltonian ( p,.10 ) can be easily diagonalized which yields the spectrum Ea^i = + dFl 



with the eigenstates 

\'^a,l)=Y.Jm-l{-^)\m) . (1.11) 

m 

Thus, all states are localized and the spectrum is the discrete Wannier-Stark ladder (|1.3| ). 

The obtained result has a transparent physical meaning. When F = Q the energy levels of 
Wannier states |/) coincide and the tunneling couples them into Bloch waves \ k) = exp(iK/)|/). 
Correspondingly, the infinite degeneracy of the level la is removed, producing the Bloch band 
eo(^) = + (Aa/2) cos(dK)J^ When F / the Wannier levels are misaligned and the tunneling 
is suppressed. As a consequence, the Wannier-Stark state involves (effectively) a finite number 
of Wannier states, as indicated by Eq. (|l.ll| ). It will be demonstrated later on that for the 
low- lying bands Eq. (|l]^) and Eq. ( 1.11| ) approximate quite well the real part of the complex 



Wannier-Stark spectrum and the resonance Wannier-Stark functions ^'^^^(x), respectively. The 
main drawback of the model, however, is its inability to predict the imaginary part of the 
spectrum (i.e. the lifetime of the Wannier-Stark states), which one has to estimate from an 
independent calculation. Usually this is done with the help of Landau-Zener theory. 

1.3. Landau-Zener tunneling 

Let us address the following question: if we take an initial state in the form of a Bloch wave 
with quasimomentum k, what will be the time evolution of this state when the external static 
field is switched on? 

The common approach to this problem is to look for the solution as the superposition of 
Houston functions 



V'(x,t) = ^c„(t)V„(x,i) , (1.12) 



1 







V'a(a:,t) = exp ( -- / dt'e„(K') ) (l)a,K'{x) , (1.13) 



where <j)a,K'{x) is the Bloch function with the quasimomentum k' evolving according to the 
classical equation of motion p = —F, i.e k' = k — Ft/h. Substituting Eq. ( pTT^ ) into the 
time-dependent Schrodinger equation with the Hamiltonian ( |1.1[) , we obtain 



= F J]X«,^(K')exp dt'[e«(K') - epi^')]^ cp , (1.14) 



^Because only the nearest ofF-diagonal elements are taken into account in Eq. (1.10), the Bloch bands are always 
approximated by a cosine dispersion relation. 
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where Xq,_^(k) = i J dxXa ^i^) d/dK,Xf3^fi{x). Neglecting the interband couphng, i.e. = 

for a 7^ /?, we have 

c^(i) for a ^ f3 and ihca = F Xa^a{n') Ca ■ (1-15) 

This solution is the essence of the so-called single-band approximation. We note that within this 
approximation one can use the Houston functions to construct the localized Wannier-Stark 



states similar to those obtained with the help of the tight-binding model. 

The correction to the solution ( |1.15| ) is obtained by using the formalism of Landau-Zener 
tunneling. In fact, when the quasimomentum k' explores the Brillouin zone, the adiabatic 
transition occurs at the points of "avoided" crossings between the adjacent Bloch bands [see, 
for example, the avoided crossing between the 4-th and 5-th bands in Fig. ^]^(a) at k = 0]. 
Semiclassically, the probability of this transition is given by 

where Aq,^^ is the energy gap between the bands and e^, e'^ stand for the slope of the bands 
at the point of avoided crossing in the limit A^^^ ^ [79|. In a first approximation, one can 
assume that the adiabatic transition occurs once for each Bloch cycle Tg = Inh/dF. Then the 
population of the a-th band decreases exponentially with the decay time 

T = h/Ta, Ta = aaFexp{-ba/F) , (1.17) 

where Oq, and ba are band-dependent constants. 

In conclusion, within the approach described above one obtains from each Bloch band a set 
of localized states with energies given by E g. (|1.3| ). However, these states have a finite lifetime 
given by Eq. ( 1.17] ). It will be shown in Sec. ^.Ij that the estimate ( |1.17| ) is, in fact, a good "first 



order" approximation for the lifetime of the metastable Wannier-Stark states. 

1.4. Experimental realizations 

We proceed with experimental realizations of the Wannier-Stark Hamiltonian (1.1 ). Originally, 



the problem was formulated for a solid state electron system with an applied external electric 
field, and in fact, the first measurements concerning the existence of the Wannier-Stark ladder 



dealt with photo- absorption in crystals |80|. Although this system seems convenient at first 
glance, it meets several difficulties because of the intrinsic multi-particle character of the system. 
Namely, the dynamics of an electron in a solid is additionally influenced by electron-phonon and 
electron-electron interactions. In addition, scattering by impurities has to be taken into account. 
In fact, for all reasonable values of the field, the Bloch time (|l]^) is longer than the relaxation 
time, and therefore neither Bloch oscillations nor Wannier-Stark ladders have been observed in 
solids yet. 



One possibility to overcome these problems is provided by semiconductor superlattices |l81| , 
which consists of alternating layers of different semiconductors, as for example, GaAs and 
AlxGai-xAs. In the most simple approach, the wave function of a carrier (electron or hole) 
in the transverse direction of the semiconductor superlattice is approximated by a plane wave 
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for a particle of mass m* (the effective mass of the electron in the conductance or valence bands, 
respectively). In the direction perpendicular to the semiconductor layers (let it be x-axis) the 
carrier "sees" a periodic sequence of potential barriers 



V{x) 



Vo 
0, 



if 3 / G Z 
else 



with 



ld\ < a/2 



(1.18) 



where the height of the barrier Vo is of the order of 100 meV and the period d ~ 100 A. 
Because the period of this potential is two orders of magnitude larger than the lattice period 
in bulk semiconductor, the Bloch time is reduced by this factor and may be smaller than the 
relaxation time. Indeed, semiconductor superlattices were the first systems where Wannier-Stark 
ladders were observed []82|-84| and Bloch oscillations have been measured in four- wave-mixing 
experiments [35,3^ as proposed in |87|. In the following years, many facets of the topics have 



been investigated. Different methods for the observation of Bloch oscillation have been applied 
1 8^ 91], and nowadays it is possible to detect Bloch oscillations at room temperature [^], to 
directly measure |^3| or even control |Q| their amplitude. Wannier-Stark ladders have been 
found in a variety of superlattice structures |9^9^, with different methods [ 100| , 101 |. The 
coupling between different Wannier-Stark ladders |102- 10^ , the influence of scattering |107- 
|109(] , the relation to the Franz-Keldysh effect 110 -112], the influence of excitonic interactions 
]11S-117] and the role of Zener tunneling ]118-121| have been investigated. Altogether, there 
is a large variety of interactions which affect the dynamics of the electrons in semiconductor 
superlattices, and it is still quite complicated to assign which effect is due to which origin. 

A second experimental realization of the Wannier-Stark Hamiltonian is provided by cold atoms 
in optical lattices. The majority of experiments with optical lattices deals with neutral alkali 



atoms such as lithium ]122], sodium ]123-125|, rubidium ]12(:-128] or cesium [129-131], but also 
optical lattices for argon have been realized ]132]. The description of the atoms in an optical 
lattice is rather simple. One approximately treats the atom as a two-state system which is 
exposed to a strongly detuned standing laser wave. Then the light-induced force on the atom is 
described by the potential [133,134] 



V{x) 



^C0s2(/CLX) 

4o 



(1.19) 



where hQji is the Rabi frequency (which is proportional to the product of the dipole matrix 
elements of the optical transition and the amplitude of the electric component of the laser field), 
kijis the wave number of the laser, and 6 is the detuning of the laser frequency from the frequency 
of the atomic transition]^ 

In addition to the optical forces, the gravitational force acts on the atoms. Therefore, a laser 
aligned in vertical direction yields the Wannier-Stark Hamiltonian 



p 

H = 1 ^ cos(2kLx) + mgx 

2m ° 



(1.20) 



where m is the mass of the atom and g the gravitational constant. An approach where one can 
additionally vary the strength of the constant force is realized by introducing a tunable frequency 

^ Th e atoms are additionally exposed to dissipative forces, which may have sub stant ial effects on the dynamics 
|l3£[ . However, since these forces are proportional to while the dipole force (1.19) is proportional to , for 
sufficiently large detuning one can reach the limit of non-dissipative optical lattices. 
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difference between the two counter-propagating waves which form the standing laser wave. If 
this difference 8lo increases hnearly in time, 8Lo{t) = 2kLat, the two laser waves gain a phase 
difference which increases quadratically in time according to 5</)(t) = kLdt^ ■ The superposition 
of both waves then yields an effective potential V{x,t) = {hVl^^/ A:5) cos^[kL{x — at^/2)], which 
in the rest frame of the potential also yields the Hamiltonian ( |1.20| ) with the gravitational force 
g substituted by a. The atom-optical system provides a much cleaner realization of the single 
particle Wannier-Stark Hamiltonian ( |1.1D than the solid state systems. No scattering by phonons 
or lattice impurities occurs. The atoms are neutral and therefore no excitonic effects have to be 
taken into account. Finally, the interaction between the atoms can be neglected in most cases 
which justifies a single particle description of the system. Indeed, Wannier-Stark ladders, Bloch 
oscillations and Zener tunneling have been measured in several experiments in optical lattices 
[1l23| , |l2l , [l29|p6| - [l3l| . 



Besides the semiconductor and optical lattices, different attempts have been made to find 
the Wannier-Stark ladder and Bloch oscillations in other systems like natural superlattices. 



optical and acoustical waveguides, etc. | 139 - 148(| . However, here we denote them mainly for 
completeness. In the applications of the theory to real systems we confine ourselves to optical 
lattices and semiconductor superlattices. 

A final remark of this section concerns the choice of the independent parameters of the systems. 
In fact, by using an appropriate scaling, four main parameters of the physical systems - the 
particle mass m, the period of the lattice d, the amplitude of the periodic potential Vq and the 
amplitude of the static force F - can be reduced to two independent parameters. In what follows 
we use the scaling which sets m = 1, Vq = 1 and d = 2tt. Then the independent parameters 
of the system are the scaled Planck constant h! (entering the momentum operator) and the 
scaled static force F' . In particular, for the system (|1.2C| ) the scaling x' = 2kLX, H' = H/Vq 
(Vb = fi'Vty4.5) gives 

i.e. the scaled Planck constant is inversely proportional to the intensity of the laser field. For 
the semiconductor superlattice, the scaled Planck constant is h' = 2T:h/ d^/mFVQ. 



1.5. This work 

In this work we describe a novel approach to the Wannier-Stark problem which has been 
developed by the authors during the last few years |149| - [164j| . By using this approach, one finds 
the complex spectrum ( |1.3D as the poles of a rigorously constructed scattering matrix. The 
suggested method is very efficient from the numerical points of view and has proven to be a 
powerful tool for an analysis of the Wannier-Stark states in different physical systems. 

The review consists of two parts. The first part, which includes chapters 2-3, deals with 
the case of a dc field. After introducing a scattering matrix for the Wannier-Stark system 
we describe the basic properties of the Wannier-Stark states, such as lifetime, localization of 
the wave function, etc., and analyze their dependence on the magnitude of the static field. A 
comparison of the theoretical predictions with some recent experimental results is also given. 
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In the second part (chapters 4-7) we study the case of combined ac-dc fields: 



H = — + V{x) + Fx + F^x cos{ujt) 
2m 



(1.22) 



We show that the scattering matrix introduced for the case of dc field can be extended to the 
latter case, provided that the period of the driving field T^^ = 2tt/uj and the Bloch period ( pTl] ) 
are commensurate, i.e. qTs = pT^ with p, q being integers. Moreover, the integer q in the 
last equation appears as the number of scattering channels. The concept of the metastable 
quasienergy Wannier-Bloch states is introduced and used to analyze the dynamical and spectral 
properties of the system ( 1.22| ). Although the method of the quasienergy Wannier-Bloch states 
is formally applicable only to the case of "rational" values of the driving frequency (in the sense 
of equation T^^/Tb = q/p)^ the obtained results can be well interpolated for arbitrary values of 

UJ. 

The last chapter of the second part of the work deals with the same Hamiltonian ( 1.22| ) 
but considers a very different topic. In chapters 2-6 the system parameters are assumed to 
be in the deep quantum region (which is actually the case realized in most experiments with 
semiconductors and optical lattices). In chapter 7, we turn to the semiclassical region of the 
parameters, where the system ( |1.22| ) exhibits chaotic scattering. We perform a statistical analysis 
of the complex (quasienergy) spectrum of the system and compare the results obtained with the 
prediction of random matrix theory for chaotic scattering. 

To conclude, it is worth to add few words about notations. Through the paper we use low 
case (f) to denote the Bloch states, which are eigenstates of the field free Hamiltonian ( |1.5D . The 
Wannier-Stark states, which solve the eigenvalue problem with Hamiltonian ( |l.lD and which are 
our main object of interest, are denoted by capital ^ . These states should not be mismatched 
with the VFannier states ( |1.8D denoted by low case ip. Besides the Bloch, Wannier, and Wannier- 
Stark states we shall introduce later on the Wannier-Bloch states. These states generalize the 
notion of Bloch states to the case of nonzero static field and are denoted by capital Thus we 
always use capital letters or <I>) to refer to the eigenfunctions for F ^ Q and low case letter 
(-0 or (j)) in the case of zero static field, as summarized in the table below. 



function 


name 




dc-field 


(t>a,K{x) 


Bloch 


delocalized eigenfunctions of the Hamiltonian Hq 


F = 




Wannier 


dual localized basis functions 


F = 




Wannier-Stark 


resonance eigenfunctions of the Hamiltonian Hw 


F/0 




Wannier-Bloch 


res. eigenfunctions of the evolution operator U{Tb) 





Chapter 2 

Scattering theory for Wannier-Stark 
systems 



In this work we reverse the traditional view in treating the two contributions of the potential to 
the Wannier-Stark Hamiltonian 

Hw = Y + ^(^) + '^(^ + 2^) = ^(^) • (2-1) 

Namely, we will now consider the external field Fx as part of the unperturbed Hamiltonian 
and the periodic potential as a perturbation, i.e. Hw = Hq + V{x), where Hq = /2 + Fx. 
The combined potential V{x) + Fx cannot support bound states, because any state can tunnel 
through a finite number of barriers and finally decay in the negative x-direction [F > 0). 
Therefore we treat this system using scattering theory. We then have two sets of eigenstates, 
namely the continuous set of scattering states, whose asymptotics define the S-matrix S{E), and 
the discrete set of metastable resonance states, whose complex energies £ = E — iT /2 are given 
by the poles of the S-matrix. Due to the periodicity of the potential V{x), the resonances are 
arranged in Wannier-Stark ladders of resonances. The existence of the Wannier-Stark ladders 



of resonances in different parameter regimes has been proven, e.g., in |25-28|. 

2.1. S-matrix and Floquet-BIoch operator 

The scattering matrix S{E) is calculated by comparing the asymptotes of the scattering states 
^s{E) with the asymptotes of the "unscattered" states ^q{E), which are the eigenstates of the 
"free" Hamiltonian 

2 

Ho = ^ + Fx, F>0. (2.2) 

In configuration space, the ^q{E) are Airy functions 

E) ~ Ai (e - 6) i-TT^O''^^ sin (C + vr/4) . (2.3) 

where ^ = ax, Co = aE/F, a = {2F/h^Y/^, and C = | (-C)^''^ Asymptotically the 

scattering states ^ s{E) behave in the same way, however, they have an additional phase shift 
ip{E), i.e. for x — > — oo we have 

^six; E) (-vr^C)-!/^ sin [ C + 7r/4 + ip{E)] . (2.4) 
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Actually, in the Stark case it is more convenient to compare the momentum space instead of the 
configuration space asymptotes. (Indeed, it can be shown that both approaches are equivalent 
|16C, 164| .) In momentum space the eigenstates ( ^.Sj ) are given by 



^o(A;;^) = exp 



6F 



Ek\ 



(2.5) 



For F > the direction of decay is the negative x-axis, so the limit k — > — oo of ^o{k; E) is the 
outgoing part and the limit A; — > oo the incoming part of the free solution. 
The scattering states ^ s[E) solve the Schrodinger equation 



Hw'i>s{E)=E^>s{E) 



(2.6) 



with Hw = Hq + V{x). (By omitting the second argument of the wave function, we stress that 
the equation holds both in the momentum and coordinate representations.) Asymptotically the 
potential V{x) can be neglected and the scattering states are eigenstates of the free Hamiltonian 
(12. 21). In other words, we have 



lim 

fe^±oo 



E) = exp 



h'^k^ Ek , , 
-6F 



(2.7) 



With the help of Eqs. (gj) and (gj) we get 



^ ' ^o{-k;E) '^s{k;E) 



which is the definition we use in the following. In terms of the phase shifts i^iE) the S-matrix 
obviously reads S{E) = exp[—i2Lp{E)] and, thus, it is unitary. 

To proceed further, we use a trick inspired by the existence of the space-time translational 



symmetry of the system, the so-called electric translation [I66|. Namely, instead of analyzing 
the spectral problem ( |2.6[ ) for the Hamiltonian, we shall analyze the spectral properties of the 
evolution operator over a Bloch period 



[/ = exp ( -^HwTb^ , Tb = ^ . 



(2. 



Using the gauge transformation, which moves the static field into the kinetic energy, the operator 
(|2.9|) can be presented in the form 



U 



(2.10) 



U 



exp , -- 



{p - Ftf 



+ V{x) 



(2.11) 



where the hat over the exponential function denotes time ordering [|. The advantage of the 
operator U over the Hamiltonian Hyi/ is that it commutes with the translational operator and, 

^Indeed, substituting into the Schrodinger equation, ihOtp/dt = Hwip, the wave function in the form tf}{x,t) = 
exp{-iFtx/h)'ip{x,t), we obtain ihdip/dt = Hwip where Hw = {p - Ftf /2 + V{x). Thus i){x,TB) = Uil'{x,0) or 
iI){x,Tb) = e'Kp{—ix)Uip{x,0). 
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thus, the formahsm of the quasimomentum can be used.0 Besides this, the evolution operator 
also allows us to treat the combined case of an ac-dc field, which will be the topic of the second 
part of this work. 

There is a one to one correspondence between the eigenfunctions of the Hamiltonian and 
the eigenfunctions of the evolution operator. Indeed, let ^s{x;E) be an eigenfunction of Hw 
corresponding to the energy E. Then the function 

^s{x; X,k) =Y^ exp(+i27r/K)^'s(x - 27r/; E) (2.12) 
I 

is a Bloch-like eigenfunction of U corresponding to the eigenvalue A = exp(— ii^Tg/Zi), i.e. 

C/$5(A, k) = A$s(A, k) , A = exp(-i^/F) . (2.13) 

Equation ( ^.131) simply follows from the continuous time evolution of the function ( 2. 12]) , which 
is ^s{x; A, K, t) = exp(+i27r/K) exp[-i(^ + 2TTFl)t/h]^ s{x - 2ttI; E), or 

$5(A, K, t) = exp(-iM/n)$5(A, K - Ft/h) . (2.14) 

Let us also note that the quasimomentum n does not enter into the eigenvalue A. Thus the 
spectrum of the evolution operator U is degenerate along the Brillouin zone. Besides this, the 
relation between energy E and A is unique only if we restrict the energy interval considered to 
the first "energy Brillouin zone", i.e. < E < 2ttF. 

When the energy is restricted by this first Brillouin zone, the transformation inverse to ( ^.12 ) 
reads 

1-1/2 

^s{E) = / dK$5(A,K) . (2.15) 

J~l/2 

This relation allows us to use the asymptotes of the Floquet-Bloch solution <I>5(A, k) instead 
of the asymptotes of the ^5'(£') in the S-matrix definition (|2.8| ). In fact, since the functions 
$5(2;; A, k) are Bloch-like solution, they can be expanded in the basis of plane waves: 

<^>s{x;\,k) =Y,Cs{n;X,K){x\n + K) , (x|n + k) = (2^)-i/2e^("+'')^ . (2.16) 

n 

From the integral ( |2.15| ) the relation {n + k|<55(A,k)) = (n + nl'^siE)) follows directly, i.e. 
in the momentum representation the functions ^s{k',E) and <I>s(A;;A,k) coincide at the points 
k = n + K. Thus we can substitute the asymptotes of ^s{k; A, k) in Eq. (|2.8[). This gives 



Q/TP\ r Cs{-n) Coin) 
n-*oo Go(— raj 05(nj 

where the energy on the right-hand side of the equation enters implicitly through the eigenvalue 
A = ex.p{—iE/F). Let us also note that by construction S{E) in Eq. ( 2.17 ) does not depend on 
the particular choice of the quasimomentum k. In numerical calculations this provides a test for 
controlling the accuracy. 



The tight-binding version of the evolution operator (2.10) was studied in Ref. [167|. 
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Figure 2.1. Matrix of the Floquet-Bloch operator U for H\y = p^/2 + cos(x) + Fx with system 
parameters h = 0.5 and F = 0.2; The absolute values of the elements are shown in a grey scale 
plot. With increasing indices the matrix tends to a diagonal one. 



2.2. S-matrix: basic equations 



Using the expansion ( ^.16|) , the eigenvalue equation ( ^.13 ) can be presented in matrix form 
Y.^!nli,nCs{n) = XCsim) , (2.18) 



where 



U^li = {m + K\U\n + k) 



(2.19) 



and the unitary operator U is given in Eq. ( 2. 111) . [Deriving Eq. ( ^.18 ) from Eq. ( |2.13| ), we took 
into account that in the plane wave basis the momentum shift operator exp(— ix) has the matrix 
elements (m| exp(— ix)|n) = 5m+i,n-] Because A does not depend on the quasimomentum k,^ 
we can set k = and shall drop this upper matrix index in what follows. For n — > ±00, the 
kinetic term of the Hamiltonian dominates the potential and the matrix U tends to a diagonal 



one. This property is exemplified in Fig. 2.1, where we depict the Floquet-Bloch matrix for the 
potential V{x) = cos(x). Suppose the effect of the off-diagonals elements can be neglected for 
Inl > A^. Then we have 



U„ 



m'-'m,n 



for Iml, Inl > 



(2.20) 



^This means that the operators exp(— ix)!/'"' are unitary equivalent - a fact, which can be directly concluded 
from the explicit form of this operator. 
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with 



(2.21) 



For the unscattered states $o('^) the formulas ( p.20 ) hold exactly for any m and, given a energy 
£^ or A = exp(— the eigenvalue equation can be solved to yield the discrete version of the 
Airy function in the momentum representation: Co{m) = exp{ih?m'^ /6F — lEm/F). With the 
help of the last equation we have 



Coin) 
Coi-n) 



exp 



2„3 



~3F~ 



2En 



(2.22) 



which can be now substituted into the S-matrix definition ( 2.17 ). 

We proceed with the scattering states $5 (A). Suppose we order the Cs with indices increasing 
from bottom to top. Then we can decompose the vector Cs into three parts, 



Cs 



( Cf^ \ 

V4-' 



(2.23) 



where C^^^ contains the coefficients for n > N, ■* contains the coefficients for n < —N — 1 

and Cg'^ contains all other coefficients for — A^— 1 < n < N. The coefficients of c|^^ recursively 
depend on the coefficient Cs{N), via 

Cs{m + 1) = (X/um+i) Cs{m) for m>N. (2.24) 

Analogously, the coefficients of ■* recursively depend on Cs{—N — 1), via 

Cs{m) = {ura+i/\)Cs{m + l) for m < -A - 1 . (2.25) 

Let us define the matrix W as the matrix [/, truncated to the size (2 A + 1) x (2 A + 1). 
Furthermore, let B]\j he the matrix W accomplished by zero column and row vectors: 



B 



N 



0* 
W 



Then the resulting equation for clP^ can be written as 



{Bn 



XI) cf 



-UN+i Cs{N + 1) e\ 



(2.26) 



(2.27) 



where is a vector of the same length as C^\ with the first element equal to one and all others 



equal to zero. For a given A, Eq. ( 2.27 ) matches the asymptotes C^^^ and ^ by linking C^\ 



via C5(A + 1) andEq. (|]2|), to and, via C5(-A^-l) and Eq. (l2^ ), to ^ Let us now 



.(-) 



introduce the row vector ei with all elements equal to zero except the last one, which equals one. 
Multiplying ei with C^'"' yields the last element of the latter one, i.e. Cs{—N — 1). Assuming 
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that A is not an eigenvalue of the matrix i?7v (this case is treated in the next section) we can 
multiply Eq. ( 2.27| ) with the inverse of {B]\f — Al), which yields 



-UN+i ei 



B 



N 



AE/F 



-1 



Csi-N - 1) 
Cs{N + 1) 

Finally, substituting Eq. ( 2.22 ) and Eq. ( |2.28| ) into Eq. ( 2.17 ), we obtain 



S{E) = lim ^(iV + l)ei 



B 



N 



-iE/F ^ 



(2.28) 



(2.29) 



with a phase factor A(N) = —uj\fCQ{N)/CQ{—N), which ensures the convergence of the limit 
N ^ oo. The derived Eq. (|2.29| ) defines the scattering matrix of the Wannier-Stark system and 
is one of our basic equations. 

To conclude this section, we note that Eq. (|2.29| ) also provides a direct method to calculate 
the so-called Wigner delay time 



TiE) 



Ah 



din S{E) 
dE 



-2h 



As shown in Ref. |152] 
h r 



t(E) = lim ^ 



d^{E) 
dE 



2(7V + 1) 



(2.30) 



(2.31) 



Thus, one can calculate the delay time from the norm of the C'g \ which is preferable to (|2.30| ) 
from the numerical point of view, because it eliminates an estimation of the derivative. In the 
subsequent sections, we shall use the Wigner delay time to analyze the complex spectrum of the 
Wannier-Stark system. 



2.3. Calculating the poles of the S-matrix 

Let us recall the S-matrix definitions for the Stark system. 



SiE) 



lim 



^s{-k;E) ^o{k;E) 



lim 

.00 ^s{k;E) ^o{-k;E) n-xx) Cs{n) Coi 

The S-Matrix is an analytic function of the (complex) energy, and we call its isolated poles 
located in the lower half of the complex plane, i.e. those which have an imaginary part less than 
zero, resonances. In terms of the asymptotes of the scattering states, resonances correspond to 
scattering states with purely outgoing asymptotes, i.e. with no incoming wave. (These are the 
so-called Siegert boundary conditions |16^ .) As one can see directly from ( 2.22| ), poles cannot 
arise from the contributions of the free solutions. In fact, Co(n)/Co(— n) decreases exponentially 
as a function of n for complex energies £ = E — ir/2. Therefore, poles can arise only from the 
scattering states C5. 

Actually, we already noted the condition for poles in the previous section. In the step from 
equation ( 2.27| ) to the S-matrix formula ( |2.29| ) we needed to invert the matrix {Bjsj — Al). We 
therefore excluded the case when A is an eigenvalue of i?7v- Let us treat it now. If A is an 
eigenvalue oi B^, the equation defining C'g^ then reads 

{Bn - Al) Cj'^ = . 



-n 



Coin) 



-n] 



(2.32) 



(2.33) 
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270 



hbar=1 ; f=0.07; kappa=0; 
NN=5; N=2*NN+1 ; jmax=16; 

dt=hbar/f/jmax; v=ones(N-1 ,1 ); 
V=0.5*(diag(v -1 )+diag(v,1 )); 
p=-hbar*([-NN:NN]'+kappa); 

U=eye(N); 
for j=1 :jmax, 
time=dt*(j-0.5); 

H=diag(0.5*(p-rtime).'^2,0)+V; 

U=expm(-i*dt*H/hbar)*U; 

end 

z=zeros(N,1 ); 

B=[z' 0; U z]; d=eig(B); 

polar(angle(d),abs(d),'*') 



Figure 2.2. The eigenvalues A of the matrix calculated for system (2A) with V{x) = cosx, 
h = 1 and F = 0.07. The numerical parameters are N = 5, jmax = 16 and k = 0. The 
eigenvalues corresponding to the first three Wannier- Stark ladders are marked by circles. On the 
right to the figure is the matlab source code which generates the depicted data. 



The scattering state Cs we get contains no incoming wave, i.e. it fulfills the Siegert boundary 
condition. In fact, the first element Cg'\N) is equal to zero, which follows directly from the 
structure of B]\[, and consequently C^^^ = 0. In addition, the eigenvalues fulfill |A| < 1,0 which 



in terms of the energy £ = E — iT/2 means F > 0. Let us also note that, according to Eq. (2.25), 
the outgoing wave Cg ^ diverges exponentially as ■'(n) ~ |A|~". 

Equation ( |2.33D provides the basis for a numerical calculation of the Wannier-Stark resonances. 
A few words should be said about the numerical algorithm. The time evolution matrix ( ^.11 ) 
can be calculated by using 2N + 1 plane wave basis states {x\n) = (27r)~-'^/^ exp(inx) via 

Jmax y . \ 

U^^^ ^Ylexpi- ^-H^^^ {tj)At ] (2.34) 
j=i 

where tj = {j- 1/2) At, At = Tb/J max and H^'^^ (tj) is the truncated matrix of the operator 
H^'^){t) = {p-Ft + hKf/2 + V{x). Then, by adding zero elements, we obtain the matrix 
Sjv and calculate its eigenvalues A. The resonance energies are given hy £ = iFlnA. As an 



example. Fig 2.2 shows the eigenvalues A^ in the polar representation for the system (|2.l|) with 



V{x) = cosx. Because of the numerical error (introduced by truncation procedure and round 



*This property follows directly from non-unitarity of Bm'. B]^Bn = 1 — e\e\. 
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error) not all eigenvalues correspond to the S-matrix poles. The "true" A can be distinguished 
from the "false" A by varying the numerical parameters N, jmax and the quasimomentum k (we 
recall that in the case of dc field A is independent of k). The true A are stable against variation 
of the parameters, but the false A are not. In Fig |2.2| , the stable A are marked by circles and can 
be shown (see next section) to correspond to Wannier-Stark ladders originating from the first 
three Bloch bands. By increasing the accuracy, more true A (corresponding to higher bands) 
can be detected. 



2.4. Resonance eigenfunctions 

According to the results of preceding section, the resonance Bloch- like functions ^a,m referred 
to in what follows as the Wannier-Bloch functions, are given (in the momentum representation) 

by 

<^aAk) =^CUn)5{n + K- k) . (2.35) 
n 

where Ca{n) are the elements of the eigenvector of Eq. ( |2.33| ) in the limit N ^ oo. The change 
of the notation $5(A,k) indicates that from now on we deal with the resonance eigen- 

functions corresponding to the discrete (complex) spectrum Sa- The Wannier-Stark states ^a,i, 
which are the resonance eigenfunction of the Wannier-Stark Hamiltonian H\y, are calculated 
by using Eq. ( 2.1^ ) and Eq. ( 2.15| ). In fact, according to Eq. ( 2.141) , the quasimomentum k of 



the Wannier-Bloch function changes linearly with time and explores the whole Brillouin zone 
during one Bloch period. Thus, one can obtain the Wannier-Stark states ^a^i by calculating the 
eigenfunction ^a,K of the evolution operator U for, say, k = and propagating it over the Bloch 
period. (Additionally, the factor exp{—i£at/h) should be compensated.) We used the discrete 
version of the continuous evolution operator, given by (|]3|) with the upper limit jmax 

substi- 



tuted by the actual number of timesteps. Resonance Wannier-Stark functions corresponding to 
two most stable resonances are shown in Fig. 



The left panel in Fig. |2.3| shows the wave functions in the momentum representation, where 
the considered interval of A; = p/h is defined by the dimension of the matrix B]\f, i.e. \k\ < N. 
The (faster than exponential) decrease in the positive direction is clearly visible. The tail in the 
negative direction reflects the decay of resonances. Although it looks to be constant in the figure, 
its magnitude actually increases exponentially (linearly in the logarithmic scale of the figure) 
as k — oo. The wave functions in the coordinate representation (right panel) are obtained 
by a Fourier transform. Similar to the momentum space the resonance wave functions decrease 
in positive x-direction and have a tail in the negative one. Obviously, a finite momentum basis 
implies a restriction to a domain in space, who's size can be estimated from energy conservation 
as \x\ < fi^N"^ /2F. Additionally the Fourier transformation introduces numerical errors due 
to which the wave functions decay only to some finite value in positive direction. We note, 
however, that for most practical purposes it is enough to know the Wannier-Stark states in the 
momentum representation. 

Now we discuss the normalization of the Wannier-Stark states. Indeed, because of the presence 
of the exponentially diverging tail, the wave functions ^oL,i{k) or ^a^i{x) can not be normal- 
ized in the usual sense. This problem is easily resolved by noting that for the non-hermitian 
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Figure 2.3. Resonance wave functions of the two most stable resonances of system (2.1) with 
parameters h = 1 and F = 0.07 in momentum and in configuration space. The ground state is 
plotted as a dashed, the first excited state as a solid line. In the second figure the first excited 
state is shifted by one space period to enhance the visibility. 



eigenfunctions (i.e. in the case considered here) the notion of scalar product is modified as 
J dxKAx)^aAx) ^ J dxf ^/x)^f/x) , (2.36) 



where ^'^^(x) and ^'^^(x) are the left and right eigenfunctions, respectively. In Fig. 2^ the 



right eigenfunctions are depicted. The left eigenfunctions can be calculated in the way described 
above, with the exception that one begins with the left eigenvalue equation c'^^ {B^ — Al) = 
for the row vector In the momentum representation, the left function ^'^^(A;) coincides 

with the right one, mirrored relative to A; = 0. (Note that in coordinate space, the absolute 
values of both states are identical.) In other words, it corresponds to a scattering state with 
zero amplitude of the outgoing wave. Since for the right wave function a decay in the positive 
A;-direction is faster than the increase of the left eigenfunction (being inverted, the same is valid 
in the negative /c-direction) , the scalar product of the left and right eigenfunctions is finite. In 
our numerical calculation we typically calculate both functions in the momentum representation 
and then normalize them according to 

dk^iAk)^^^Ak) = {^a,l\^,3,n) = S'^p . (2.37) 

(Here and below we use the Dirac notation for the left and right wave functions.) Let us also 
recall the relations 

^a,iix) = *a,o(x - 27r0 (2.38) 
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Figure 2.4. Comparison of the wave functions calculated within the different approaches for 
h = 2 and F = 0.01,0.03,0.1, shown on a linear (top) and on a logarithmic scale (bottom). The 
dotted line is the tight-binding, the dashed line the single-band and the solid line is the scattering 
result. 



for the wave functions in the coordinate representation and 

= exp(i27rZA;)^'«,o(^) (2.39) 

in the momentum space. Thus it is enough to normaUze the function for / = 0. Then the 
normahzation of the other functions for / ^ wiU hold automatically. For the purpose of future 
reference we also display a general (not restricted to the first energy Brillouin zone) relation 
between the Wannier-Bloch and Wannier-Stark states 



/■1/2 

^a,i= dKexp(-i27r/K)$„,« (2.40) 

J -1/2 



(compare with Eq. (|1.8|). 
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It is interesting to compare the resonance Wannier-Stark states with those predicted by the 
tight-binding and single-band models. Such a comparison is given in Fig. p.4| , where the ground 
Wannier-Stark state for the potential V{x) = cos x is depicted for three different values of the 
static force F. As expected, for small F, where the resonance is long-lived, both approximations 
yield a good correspondence with the exact calculation. (In the limit of very small F the 
single-band model typically gives a better approximation than the tight-binding model.) In the 
unstable case, where the resonance state has a visible tail due to the decay, the results differ in 
the negative direction. On logarithmic scale one can see that the order of magnitude up to which 
the results coincide is given by the decay tail of the resonances. In the positive x-direction the 
resonance wave functions tend to be stronger localized. It should be noted that in Fig. 2.4 we 



considered the ground Wannier-Stark states only for moderate values of the static force F < 0.1. 
For larger F, because of the exponential divergence, the comparison of the resonance Wannier- 
Stark states with the localized states of the single-band model loses its sense. The same is also 
true for higher (a > 0) states. Moreover, the value of F, below which the comparison is possible, 
rapidly decreases with increase of band index a. 



Chapter 3 

Interaction of Wannier-Stark ladders 



In this chapter we give a complete description of the dependence of the width F of the Wannier- 
Stark resonances on the parameters of the Wannier-Stark Hamiltonian. In scaled units, the 
Hamiltonian has two independent parameters, the scaled Planck constant h and the field strength 
F. In our analysis we fix the value of h and investigate the width as a function of the field 
strength. The calculated lifetimes r = h/T are compared with the experimentally measured 
lifetimes of the Wannier-Stark states. 

3.1. Resonant tunneling 

To get a first glimpse on the subject, we calculate the resonances for the Hamiltonian ( ^.1| ) 
with V{x) = cosx for h = 1. For the chosen periodic potential the field-free Hamiltonian has 
two bands with energies well below the potential barrier. For the third band, the energy e2{K) 
can be larger than the potential height. Therefore, with the field switched on, one expects two 
long-lived resonance states in each potential well, which are related to the first two bands. 

Figure |3.1| (a) shows the calculated widths of the six most stable resonances as a function 
of the inverse field strength 1/F. The two most stable resonances are clearly separated from 
the other ones. The second excited resonance can still be distinguished from the others, the 
lifetime of which is similar. Looking at the lifetime of the most stable state, the most striking 
phenomenon is the existence of very sharp resonance-like structures, where within a small range 
of F the lifetime can decrease up to six orders of magnitude. In Fig. |3.1| (b), we additionally 
depict the energies of the three most stable resonances as a function of the inverse field strength. 
As the Wannier-Stark resonances are arranged in a ladder with spacing AE = 2itF, we show 
only the first energy Brillouin zone < E/F < 2tt. Let us note that the mean slope of the 
lines in Fig. |3.1| (b) defines the absolute position E^ of the Wannier-Stark resonances in the limit 
F — > 0. As follows from the single band model, these absolute positions can be approximated 
by the mean energies of the Bloch bands. Depending on the value of E^, we can identify a 
particular Wannier-Stark resonance either as under- or above-barrier resonance.^] 

Comparing Fig. |3.1| (b) with Fig. |3.1| (a), we observe that the decrease in lifetime coincides 
with crossings of the energies of the Wannier-Stark resonances. All three possible crossings 

^This classification holds only in the limit F ^ 0. In the opposite limit all resonances are obviously above-barrier 
resonances. 
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Figure 3.1. a) Resonance width of the 6 most stable resonances as a function of the inverse 
field strength 1/F. b) Energies of the 3 most stable resonances as a function ofl/F (solid line: 
most stable resonance, dashed line: first excited resonance, dashed dotted line: second excited 
resonance). Parameters are V{x) = cosx and h= 1. 
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Figure 3.2. Wannier-Stark resonances in different minima of the potential V{x) = cos(2;) +Fx: 
The most stable resonance and some members of the first excited Wannier-Stark ladder are 
shown. The parameters are H = 1.0 and F = 0.08. 



manifest themselves in the hfetime: Crossings of the two most stable resonances coincide with 
the sharpest peaks in the ground state width. The smaller peaks can be found at crossings of the 
ground state and the second excited state. Finally, crossings of the first and the second excited 
state fit to the peaks in the width of the first excited state. The explanation of this effect is the 
following: Suppose we have a set of resonances which localize in one of the 27r-periodic minima 
of the potential V{x) = cosx + Fx. Let AE'^^j = Ea — be the energy difference between 
two of these states. Now, due to the periodicity of the cosine, each resonance is a member of a 
Wannier-Stark ladder of resonances, i.e. of a set of resonances with the same width, but with 
energies separated by Ai? = 2'kF. Figure 3.2 shows an example: The two most stable resonances 
for one potential minimum are depicted, furthermore two other members of the Wannier-Stark 
ladder of the first excited resonance. To decay, the ground state has to tunnel three barriers. 
Clearly, if there is a resonance with nearly the same energy in one of the adjacent minima, 
this will enhance the decay due to phenomenon of resonant tunneling. The strongest effect will 
be given for degenerate energies, i.e. for liiFl = AEa^/3, which can be achieved by properly 
adjusting F, because the splitting AEa,/3 ^ E^ — £'* is nearly independent of the field strength. 
For the case shown in Fig. 3.2, such a degeneracy will occur, e.g., for a slightly smaller value 



F ~ 1/14.9 (see Fig. 3J). Then we have two resonances with the same energies, which are 



separated by two potential barriers. In the next section we formalize this intuitive picture by 
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introducing a simple two-ladder model. 



3.2. Two interacting Wannier-Stark ladders 

It is well known that the interaction between two resonances can be well modeled by a two- 



state system [34,16E-171|. In this approach the problem reduces to the diagonalization of a 2 x 2 
matrix, where the diagonal matrix elements correspond to the non-interacting resonances. In 
our case, however, we have ladders of resonances. This fact can be properly taken into account 
by introducing the diagonal matrix in the form [155,160| 



Uo = exp ( -1 — 



En 



- iro/2 
El 





-iri/2 



(3.1) 



It is easy to see that the eigenvalues Ao,i(-F) = exp[— i(ii^o,i — iro,i/2]/i^) of Uq correspond to 
the relative energies of the Wannier-Stark levels and, thus, the matrix Uq models two crossing 
ladders of resonances .0 Multiplying the matrix Uq by the matrix 



int 



exp 



le 



1 

1 



cos e 1 sm e 
i sin e cos e 



(3.2) 



we introduce an interaction between the ladders. The matrix UoUint can be diagonalized ana- 
lytically, which yields 



, Ao + Ai 

A± = cos e ± 



Ao + Ai 



cos e 



AqAi 



1/2 



X± = exp — i 



E± - ir±/2 

F 



(3.3) 



Based on Eq. ( |3.3|) we distinguish the cases of weak, moderate or strong ladder interaction. 

The value e = obviously corresponds to non-interacting ladders. By choosing e ^ but 
e <^ 7r/2 we model the case of weakly interacting ladders. In this case the ladders show true 
crossing of the real parts and "anticrossing" of the imaginary parts. Thus the interaction affects 
only the stability of the ladders. Indeed, for e <C 7r/2 Eq. ( |3.3D takes the form 



Ah 



A, 



0,1 



1 ± 



e2 AojfAi 
2 Ai — Ao 



(3.4) 



It follows from the last equation that at the points of crossing (where the phases of Aq and Ai 
coincide) the more stable ladder (let it be the ladder with index 0, i.e. Fq < Fi or |Ao| > |Ai|) 
is destabilized (|A+| < |Ao|) and, vice versa, the less stable ladder becomes more stable (|A_| > 



I All). The case of weakly interacting ladders is illustrated by the left column in Fig. 3.3 
By increasing e above 



sm £(. 



|Ao| -|Ai| 
|Ao| + |Ai| 



(3.5) 



the case of moderate interaction, where the true crossing of the real parts E± is substituted 
by an anticrossing, is met. As a consequence, the interacting Wannier-Stark ladders exchange 

^ The resonance energies in Eq. ( [3.l| ) actually depend on F but, considering a narrow interval of F, this dependence 
can be neglected. 
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Figure 3.3. Illustration to the two-ladder model. Parameters are £q = 0.3 — il.l • 10~^, £i = 
0.8 — iO.9 • 10~^, and e = 0.2 (left column), e = 0.4 (center), and e = ■k/2 — 0.1 (right column). 
Upper panels show the energies E±, lower panels the widths V±. 



their stability index at the point of the avoided crossing (see center column in Fig. |3.3| ). The 
maximally possible interaction is achieved by choosing e = 7r/2. Then the eigenvalues of the 
matrix C/o^int are \± = ibi(AoAi)^/-^ which corresponds to the "locked" ladders 

E± = {Eo + Ei)/2 ± 7rF/2 , r± = (Fq + Fi)/2 . (3.6) 

In other words, the energy levels of one Wannier-Stark ladder are located exactly in the middle 
between the levels of the other ladder (right column in Fig. |3.3| ). 

3.3. Wannier-Stark ladders in optical lattices 

In the following two sections we give a comparative analysis of the ladder interaction in optical 
and semiconductor superlattices. It will be shown that the character of the interaction can be 
qualitatively deduced from the Bloch spectrum of the system. 

We begin with the optical lattice, which realizes the case of a cosine potential (see Sec. ^^). 
A characteristic feature of the cosine potential is an exponential decrease of the band gaps as 



-E — > oo [see Fig. 1.2(a), for example]. In order to get a satisfactory description of the ladder 



interaction for F 7^ 0, it is sufficient to consider only the under-barrier resonances and one or 
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two above-barrier resonances. In particular, for the parameters of Fig. 3.1 it is enough to "keep 



track" of the resonances belonging to the first three Wannier-Stark ladders. It is also seen in 



Fig. 3.1 that the case of true crossings of the resonances is realized almost exclusively, i.e. the 
ladders are weakly interacting (which is another characteristic property of the cosine potential). 
The behavior of the resonance widths (F) at the vicinity of a particular crossing is captured 
by Eq. (p.4|). Moreover, extending the two-ladder model of the previous section to the three 
ladder case and assuming the coupling constants in the form 

€a = aaexp{-ba/F) , (3.7) 

(which is suggested by the semiclassical arguments of Sec. p!^ ) the overall behavior of the 
resonance width can be perfectly reproduced (see Fig. |3.^). The procedure of adjustment of the 



model parameters Oq, and ba is carefully described in Ref. |160| . 

The lifetime of the Wannier-Stark states (given by r = h/Ta) as the function of static force 



was measured in an experiment with cold sodium atoms in a laser field |124]. The setting of the 
experiment [[L24| yields the accelerated cosine potential (the inertial force takes the role of the 
static field) and an effective Planck constant h = 1.671. For this value of the Planck constant 



one has only one under-barrier resonance, and the two-ladder model of Sec. |3.2| is already a 



good approximation of the real situation. Figure 3.5 compares the experimental results for the 
lifetime of the ground Wannier-Stark states with the theoretical results. The axes are adjusted 
to the experimental parameters. Namely, the field strength in our description is related to the 
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acceleration (km/s ) 



Figure 3.5. Lifetime of the ground Wannier state as a function of the external field. The solid 



line is the theoretical prediction, the circles are the experimental data of Ref. , 124 1- The insert 
blows up the interval 4000m/s^ < a < lOOOOm/s^ considered in the cited experiment. 



acceleration in the experiment by the formula F ~ 0.0383a, where a is measured in km/s'^, 
and the unit of time in our description is approximately 1.34/xs. The experimental data follow 
closely the theoretical curve. (Explicitly, the analytical form of the displayed dependence is 
given by Eq. ( ^^ ) with e = aexp(— a = 1.0, b = 0.254.) In particular we note that the 
theory predicts a local minimum of the lifetime at a = 5000m/s^, which corresponds to the 
crossing of the ground and the first excited Wannier levels in neighboring wells. Unfortunately, 
the experimental data do not extend to smaller accelerations, where the theory predicts much 
stronger oscillations of the lifetime. 

3.4. Wannier-Stark ladders in semiconductor superlattices 



We proceed with the semiconductor superlattices. As mentioned in Sec. 1.4, the semiconductor 



superlattices are often modeled by the square-box potential ( 1.18 ), where a and b = d — a are 
the thickness of the alternating semiconductor layers. For the square-box potential ( |1.18| ) the 
width of the band gaps decreases only inversely proportional to the gap's number. Because of 
this, one is forced to deal with infinite number of interacting Wannier-Stark ladders. However, 



as was argued in Ref. [163|, this is actually an over-complication of the real situation. Indeed, 



the potential ( 1.18| ) is only a first approximation for the superlattice potential, which should be 



a smooth function of x. This fact can be taken into account by smoothing the rectangular step 
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in ( p^ as 

V{x) = tanh(T(x + a-K/2d) - tanh(T(x - a-K/2d) - 1 (3.8) 

for example. (Here we use scaled variables, where the potential is 27r-periodic and |V^(x)| < 1.) 
The parameter cr"^ defines the size of the transition region between the semiconductor layers 
and, in natural units, it cannot be smaller than the atomic distance. The smoothing introduces 
a cut-off in the energy, above which the gaps between the Bloch bands decrease exponentially. 
Thus, instead of an infinite number of ladders associated with the above-barrier resonances, we 
may consider a finite number of them. The interaction of a large number of ladders originating 



from the high-energy Bloch bands was studied in some details in Ref. |163]. It was found that 
they typically form pairs of locked [in the sense of Eq. ( |3.6| )] ladders which show anticrossings 
with each other. 

Since the lifetime of the above-barrier resonances is much shorter than the lifetime of the 
under-barrier resonances one might imagine that the former are of minor physical importance. 
Although this is partially true, the above-barrier resonances cannot be ignored because they 
strongly affect the lifetime of the long-lived under-barrier resonances. This is illustrated in 
Fig. |3.6| , where the resonance structure of the Wannier-Stark Hamiltonian with a periodic po- 
tential given by Eq. ( |3.8| ) and h = 3.28 is depicted as a gray-scaled map of the Wigner delay 
time (|]3|). In terms of Fig. U, this way of presentation of the numerical results means that 
each line in the lower panel has a "finite width" defined by the value F in the upper panel. In 



fact, asssuming a Wigner relation [199| we get 

/ oo 

r{E)=To + Y.[ Im 



£a + 2ttFI - E 



(3.9) 



where each term in the sum over a is just a periodic sequence of Lorentzians with width Fq,. 
(We recall that, by definition, t{E) is a periodic function of the energy.^ In the case of a 
large number of interacting ladders (i.e. in the case currently considered here, where more 
than 10 above-barrier resonances contribute to the sum over a) we find this presentation more 
convenient because it reveals only narrow resonances, while the wide resonances contribute to 
the background compensated by the constant tq. For the chosen value of the scaled Planck 
constant, h = 3.28, the periodic potential ( ^ ) supports only one under-barrier resonance, seen 
in the figure as a broken line going from the upper-left to the lower-right corners. Wide above- 
barrier resonances originating from the second and third Bloch bands and showing anticrossings 
with the ground resonance can be still identified, but the other resonances are indistinguishable 
because of their large widths. Nevertheless, the existence of these resonances is confirmed 
indirectly by the complicated structure of the "visible lines" . 

In conclusion, in comparison with the optical lattices, the structure of the Wannier-Stark 
resonances in semiconductor superlattices is complicated by the presence of large number of 
above-barrier resonances. Besides this, in the semiconductor superlattices a strong interaction 
between the ladders is the rule, while the case of weakly interacting ladders is typical for optical 
lattices. 



^The quantity (3.9) can be also interpreted as the fluctuating part of tfie (normalized) density of states of the 
system. 
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Figure 3.6. Grey-scaled map of the Wigner delay time ^3. f\ ) for the smoothed square-box potential 
(3.6 ). The parameters are h = 3.28, a/b = 39/76 and 5 = 0.25. 



Chapter 4 

Spectroscopy of Wannier-Stark 
ladders 



In this chapter we discuss the spectroscopy of Wannier-Stark ladders in optical and semiconduc- 
tor superlattices. We show how the different spectroscopic quantities (measured in a laboratory 
experiment) can be directly calculated by using the formalism of the resonance Wannier-Stark 
states. 

4.1. Decay spectrum and Fermi's golden rule 

The spectroscopy approach assumes that one probes a quantum system by a weak ac field 
Fi^x cos{ujt) with tunable frequency uj. In our case, the system consists of different Wannier- 



Stark ladders of resonances, the two most stable of which are schematically depicted in Fig. |4J. 
The driving induces transitions between the ground and the excited statesQ. Scanning the 
frequency lo sequentially activates the different transition paths and the different Wannier states 
of the excited ladder are populated. Because the excited states are typically short-lived, they 
decay before the driving can transfer the population back to the ground state, i.e. before a Rabi 
oscillation is performed. Then the decay rate of the ground state is determined by the transition 
rate D{uj) to the excited Wannier-Stark ladder. The width is written as 

roM=ro + DM, (4.1) 

where Fq takes into account the decay in the absence of driving. In what follows we shall refer 
to the quantity ro(a;) as the induced decay rate or the decay spectrum. In Sec. ^ we calculate 
the induced decay rate rigorously by using the formalism of quasienergy Wannier-Stark states. 
It will be shown that the decay spectrum is given by 



FoM=ro + ^^^Im 



2 

/3>0 L 



{Ep^i + 27rFL - Eo^i - hw) - 



(4.2) 



^Actually, transitions within the same ladder are also induced, but their effect is important only for u ^ ujb ~ 
2-KF/h. Here we shaU mainly consider the case a; S> (j-'s, where the transitions within the same ladder can be 
ignored. 
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Figure 4.1. Schematic illustration of the transitions induced by a periodic driving. The positions 
of the ground and the first excited Wannier-Stark ladder are shown for F = 0.04 and h = 1.5. 
The width of the states is symbolized by the different strength of the lines. 



where and uj are the amplitude and frequency of the probing field and 
VqA^) = (^0,zk|^'/3,/+L)(^'/3,i+L|x|^'o,/) 



(4.3) 



is the square of the dipole matrix element between an arbitrary ground Wannier-Stark state 
^'o,/(a;) and the upper Wannier-Stark state ^fs^i+ii^) shifted by L lattice periods. We would 
like to stress that, because for the resonance wave functions (^'q,^/|x|*I'/3^/') ^ {"^ f3^ii\x\'^a,i)* , the 
square of the dipole matrix element Iq^^(-L) is generally a complex number. 

it is useful to discuss its relation to Fermi's 



(4.4) 



To understand the physical meaning of Eq. (|4.2| 
golden rule, which reads 



dE 



dx^*E{^)x^Eo, (x) p{E)6{E - Eo,i - hw) 



in the notations used. In Eq. ( |4.4| ), the ^six) are the hermitian eigenfunctions of the Hamil- 
tonian (2^) (i.e., E is real and continuous) and p{E) is the density of states. For the sake of 
simplicity we also approximate the ground Wannier-Stark resonance by the discrete level Eqi. 



Then Eq. (4^) describes the decay of a discrete level into the continuum. Assuming, for a mo- 
ment, that the continuum is dominated by the first excited Wannier-Stark ladder, the density 
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of states p{E) is given by a periodic sequence of Lorentzians with width Ti, i.e. 



Substituting the last equation into Eq. ( |4.4| ) and integrating over E we have 



(4.6) 



In the case Fi <C 27rF the Lorentzians in the right-hand side of Eq. ( [4.6| ) are (5-hke functions of 
the argument fuj = Ei^i + 2'kFL — Eq^i. Thus the transition matrix element can be moved under 
the summation sign, which gives 



where 



|yo,/3p(i) = I y" dx^%^^^^^{x)x^E,A^)^ (4.8) 

(here we again included the possibility of transitions to the higher Wannier ladders, which is 
indicated by the sum over /3). It is seen that that the obtained result coincides with Eq. (|4.2| ) if 
the coefficients | Vo,/3p(-^) ^'^^ identified with the squared dipole matrix elements ( [4.31 ). Obviously, 
this holds in the limit F — > 0, when the resonance wave functions can be approximated by the 
localized states. For a strong field, however, Eq. ( |4.7| ) is a rather poor approximation of the 
decay spectrum. In particular, it is unable to predict the non-Lorentzian shape of the lines, 
which is observed in the laboratory and numerical experiments and which is correctly captured 
in Eq. ( |4.2|) by the complex phase of the squared dipole matrix elements VqJL). 

To proceed further, we have to calculate the squared matrix elements ( [4. 3D . A rough esti- 
mate for Vqp^L) can be obtained on the basis of Eq. (|1.11| ), which approximates the resonance 
Wannier-Stark state by the sum of the localized Wannier states: ^a,i = Jm-i{^a/'i'^F)'ipa,m- 



The typical experimental settings (see Sec. [4. 3]) correspond to Aq/AttF <C 1 and Aj^/ATrF > 1. 



Then the values of the matrix elements are approximately 

Vi^piL) ^ \Vo,M\L) ^ \{^|;o,l\x\^p,l)\'Jl (^) , (4.9) 

which contribute mainly in the region L < Ap/AirF, the localization length of the excited 
Wannier-Stark states. The degree of validity of this result is discussed in the next subsection. 

4.2. Dipole matrix elements 

In this subsection we calculate the dipole matrix elements 

VaA^ - I') = {^aM^p,l') (4.10) 
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beyond the tight-binding approximation. We shall use Eq. ( |2.4[l| ) 

^a,iix) = j dKe-'2^'''$«,«(x) , ^aM) = e'^^'XaAx) , XaA^) = XaA^ + 27r) , (4.11) 

which relates the Wannier-Stark states '^a,i{x) to the Wannier-Bloch states ^a^^i^)- As follows 
from the results of Sec. |2|, the function Xo,k(^) generated from Xafi{x) by propagating it 

in time 

IXa,K) = exp {\ ^) C/(t)|Xa,o) , (4.12) 
where U{t) is the continuous version of the operator U defined in Eq. ( |2.11| ) and the quasi- 



momentum k is related to time t by k = —Ft/h. Substituting Eq. ( 4.11| ) and Eq. ( 4.12 ) into 
Eq. (4.10) we obtain the dipole matrix elements as the Fourier image 



y„,^(Z - I') = 27TI 6'^''^ + I dKe^2-('-'>x,,^(/t) (4.13) 

of the periodic function 

d 1 (p + hnf' £a 

Xa,p{n) = i {Xa,K\ -Q^Xf3,K) = -p {XaA ^ h V {x) |x/3,k) - y- 5a,f3 ■ (4.14) 

The last two equations provide the basis for numerical calculation of the transition matrix 
elements. We also recall that one actually needs the square of the matrix elements ( [4. 3] ) but 
not the matrix elements themselves (which are defined up to an arbitrary phase). Thus we first 
calculate V^^^(L) and „(L) for L = 0, ±1, . . . and then multiply them term by term. 

In Fig. |4.2| we depict the squared dipole matrix elements between the ground and first excited 
Wannier-Stark states for V{x) = cos x, a moderate values of the static force F = 0.04 and values 
of the scaled Planck constant in the interval 1 < h < 2.5. For h = 1.0 the Bloch bands width 
Ai ~ 0.05 is much smaller than 4:TtF ~ 0.5 and the upper Wannier-Stark state is essentially 
localized within single potential well.^ Then only "vertical" transitions, L = 0, are possible 
between the ground and first excited Wannier ladders. By increasing h the localization length 
of the upper state grows (proportional to the band width) and more than one matrix element 
may differ from zero. Simultaneously, the Wannier levels move towards the top of the potential 
barrier (for ^ > 1.6 the upper Wannier level is already above the potential barrier) and the 
Wannier state looses its stability (Fi = 1.90 • 10"^^, 1.35 • lO'^, 5.24 • lO'^, and 1.14 • 10'^, for 
h = 1, 1.5, 2, and 2.5). Because for short-lived resonances the tight-binding result ( |1.11| ) is a 
rather poor approximation of the resonance wave functions, we observe an essential deviation 
from Eq. ([4.g|). In particular we note a strong asymmetry of the matrix elements with respect 
to L. It appears that the transitions "down the ladder" are enhanced in comparison with the 
transitions "up the ladder" . At the same time, for weak far transitions (L ^ 1) the situation is 
reversed [see Fig. B^d) and Fig. |4l(b) below]. 



Substituting the calculated matrix elements into Eq. (4^), we find the decay spectra of the 



system. The solid line in Fig. 4.3 shows the decay spectra for h = 1.5, 2.0, 2.5. As expected, 



ro(c<-') has number of peaks with the same width Fi separated by the Bloch frequency ujb- The 
^The ground Wannier-Stark state is localized within one well for all considered values of the scaled Planck constant. 
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Figure 4.2. The absolute values of the squared dipole matrix elements (^) for V{x) = cosx, 
F = 0.04 andh=l (a), 1.5 (b), 2 (c), and 2.5 (d). 



relative heights of the peaks are obviously given by the absolute values of the squared dipole 
matrix elements shown in Fig. |4.2| , while the shape of the lines is defined by the phase of V^^^ (L) . 
As mentioned above, the phases of the squared dipole matrix elements are generally not zero 
and, therefore, the shape of the lines is generally non-Lorentzian. In other words, we meet 
the case of Fano-like resonances |172|. For the sake of comparison the dashed lines in Fig. 4.3 



show the results of an exact numerical calculation of the decay rate. A good correspondence 
is noticed. The discrepancy in the region of small driving frequency is due to the rotating 
wave approximation (which is implicitly assumed in the Fermi golden rule) and the effect of the 
diagonal matrix elements V^^{L) (which are also ignored in the Fermi golden rule approach). 
In principle, the region of small driving frequency requires a separate analysis. 

In conclusion, we discuss the effect of direct transitions to the second excited Wannier ladder. 
For the case h = 2 the squared dipole matrix elements Vqi{L) and VQ2iL) are compared in the 
left column of Fig. [4.4| . It is seen that the main lines in Fig. [4.4| (c) are ten times smaller than 



those in Fig. 4.4(a). Thus the effect of higher transitions can be neglected. We note, however, 
that this is not always the case. In the next section we consider a situation when the direct 
transitions to the second excited Wannier ladder can not be ignored. 
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Figure 4.3. Comparison of the lj- dependence in ( [^.^ j (solid line) with the exact numerical 
calculation of the induced decay rate (dashed line). Parameters are F = 0.04, F^^ = 0.02 and 
h = 1.5 (left), h = 2.0 (middle) and h = 2.5 (right panel). 



4.3. Decay spectra for atoms in optical lattices 

The induced decay rate ro(a;) was measured for the system of cold atoms in the accelerated 



standing laser wave [|123| , |125| . Because the atoms are neutral, the periodic driving of the system 



was realized by means of a phase modulation of the periodic potential: 

if = y + cos[x + ecos(wt)] +Fx , (4.15) 



Using the Kramers-Henneberger transformation [ 173 - 176 | ^ the Hamiltonian (|4.15| ) can be pre- 
sented in the form 

H = Y + cos(x) + Fx + F^x cos{Lot) , F^ = ecu^ . (4.16) 

Thus, the phase modulation is equivalent to the effect of an ac field. Considering the limit of 
small e, where cos[x + e cos{Lut)] cosx + esinxcos(tiJt), we can adopt Eq. ( [4. 2D of the previous 
section to cover the the case of phase modulation. Namely, the amplitude in Eq. (^) should 
be substituted by e and the squared dipole matrix elements ( |4.3| ) by the squared matrix elements 

Wlp{L) = {^o^i\ sin x\^f3,i+L){'^f3,i+L\ sinxl^'o,;) . (4.17) 



^The Kramers-Henneberger transformation is a canonical transformation to the oscillating frame. In the classical 
case it is defined by the generating function J-{p', x, t) — [p' + eui sm{ujt)] [x + e cos(ajt)] . In the quantum case one 
uses a substitution %j}{x,t) = exp[— j_FL sm{uji)x /huj)il}{x,t) together with the transformation x' = x + ecos(cjt). 
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Figure 4.4. The absolute values of the transition coefficients Vqi{L) (a), Wqi{L) (h), Vq2{L) 
(c), and W^^iL) (d). 



Moreover, according to the commutator relation for the Hamiltonian of the non-driven system 
h-'^[HwAHw,x\] = -s\nx + F , (4.18) 
the squared matrix elements Wq ^{L) are related to the squared dipole matrix elements VQp{L) 

by 



0,1 



h 



(4.19) 



It follows from the last equation that the way of driving realized in the optical lattices suppresses 
the transition down the ladder and enhances the transition up the ladder. This is illustrated 
in Fig. 4.4, where we compare the squared matrix elements VFq^(L) and VJ)^^(L) for /3 = 1,2 
calculated on the basis of Eq. ( 4.17] ) and Eq. ( [l.sp , respectively. It is seen that the practically 
invisible tail of far transitions in Fig. ^^(a) shows up in Fig. [4. 4K b). Besides this, for L ^ 1 the 
squared matrix elements between the ground and second excited Wannier-Stark states are larger 
than those between the ground and first excited one. Because the width of the second excited 
Wannier-Stark resonance r2 is lager than Fi (and actually larger than the Bloch energy), the 
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Figure 4.5. Decay spectra as a function of the driving frequency oj. Parameters are F = 0.04, 
£ = 0.02 and h = 1.5 (left), h = 2.0 (middle) h = 2.5 (right panel). The exact numerical 
calculation (dashed lines) are compared to the model prediction (solid lines). Note a complicated 
structure of the decay spectra in the high-frequency region caused by the interference of the 
transitions to the first and second excited Wannier-Stark ladders. 



transition to the first and second excited Wannier ladders may interfere. Indeed, tliis is the case 



usuahy observed in the high-frequency regime of driving (see Fig. 4.5, which should be compared 
with Fig. HJ). 

We proceed with the experimental data for the spectroscopy of atomic Wannier-Stark ladders 
|123| ] (note also the improved experiment | 125 |). The setup in the experiment [123| is as follows. 
Sodium atoms were cooled and trapped in a far-detuned optical lattice. Then, introducing a 
time-dependent phase difference between the two laser beams forming the lattice, the lattice 



was accelerated (see Sec. L4). After some time, only atoms in the ground Wannier-Stark states 
survived, i.e. a superposition of ground ladder Wannier-Stark states was prepared. Then an 
additional phase driving of frequency uj was switched on and the survival probability. 



Pt{uj) = exp 



Toiu;)t 

h 



(4.20) 



was measured. The experiment was repeated for different values of uj. In scaled units the 
experimental settings with Vo/h = 75 it 7kHz (we choose the value Vo/h = 68kHz, which 
is used in all numerical simulations in [ |123| ) and a = 1570m/s^ correspond to h = 1.709 
and F = 0.0628. (For these parameters the ground and first excited state have the widths 
Fq = 2.38 • 10~^ and Fi = 6.11 • 10~^, respectively.) The timescale in the experiments is 1.37//S, 
and the Bloch frequency is u^bI'^t^ = 26.85kHz. The driving amplitude was e = 0.096. The 



left panel of Fig. 4.6 shows the decay spectra as a function of the frequency in this case. The 
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Figure 4.6. The left panel shows the induced decay rate of the ground Wannier-Stark state as 
function of the driving frequency for h = 1.709, F = 0.0628 and e = 0.096. The right panel 



compares the experimental data from [125 J with the calculated survival probability Pt{uj) for 
t = 300/xs. 



vertical transition dominates the figure, accompanied by the two transitions with L = ±1 and 
a tail of transitions with positive L ^ 1. In the right panel, the experimental data for the 
survival probability Pt{uj) are compared to our numerical data. The time t is taken as an 
adjustable parameter and chosen such that the depth of the peaks approximately coincide. The 
curve shows the survival probability at t = 300//S corresponding to t = 219 in scaled units. A 
good correspondence between experiment and theory is noticed. The minima of the survival 
probability appear when the driving frequency fits to a transition. The relative depth of the 
minima reflecting the size of the transition matrix elements agrees reasonably. Furthermore, 
the asymmetric shape of the minimum between Aojb and 5ujb is reproduced. Note that the 
experimental data also allow to extract the width of the first excited state from the width of the 
central minimum: Fi 0.3wB 6.9-10-2, which is in reasonable agreement with the numerical 
result Fi = 6.11 • lO^^. 
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4.4. Absorption spectra of semiconductor superlattices 



Equation (|4.2|) of Sec. 4.1 can be generalized to describe the absorption spectrum D{u>) of 



undoped semiconductor superlattices [163|. This generalization has the form 



a,f3 L 



- Ki + ^dFL + Eg-nw)- i{T% + r^)/2 



(4.21) 



where the upper indices e and h refer to the electron and hole Wannier-Stark states, respectively, 
Eg is the energy gap between the conductance and valence bands in the bulk semiconductor, 
and 



(4.22) 



is the square of the overlap integral between the hole and electron wave functions. Repeating 
the arguments of Sec. iA it is easy to show that in the low- field limit Eq. ( 4.21J ) is essentially 
the same as the Fermi golden rule equation 



D{u:) 



j j dE^dE^\ j dx^'{x;E^)^^{x;E^)\%''{E^)p''{E^)6{E^ -E^ + Eg-hLj) , (4.23) 



where p^{E'^) and p^{E^) are the one-dimensional electron and hole densities of states. According 
to Ref. i7|,|T2l the quantity D{uj), which can be interpreted as the probability of creating the 
electron- hole pair by a photon of energy huj (the electron-hole Coulomb interaction is neglected) , 
is directly related to the absorption spectrum of the semiconductor superlattices measured in 
the laboratory experiments. 

It follows from Eq. ( [4.21| ) that the structure of the absorption spectrum depends on the values 
of the squared overlap integral Eq. ( [4.22 ) which, in turn, depend on the value of the static field. 
In the low-field regime the Wannier-Stark states are delocalized over several superlattice periods 
and many transition coefficients jsiL/) differ from zero. In the high- field regime the Wannier- 
Stark states tend to be localized within a single well and the vertical transitions L = become 
dominant. We would like to stress, however, that the process of localization of the Wannier- 
Stark states is always accompanied by a loss of their stability. As mentioned above, the latter 
process restricts the validity of the tight-binding results concerning a complete localization of 
the Wannier-Stark states in the limit of strong static field. 

As an illustration to Eq. (4.22), Fig. 4.7 shows the absorption spectrum of the semiconductor 
superlattice studied in the experiment [121|.0 (This should be compared with the absorption 
spectrum calculated in Ref. |47| by using a kind of finite-box quantization method.) The depicted 
result is a typical example of a Wannier-Stark fan diagram. By close inspection of the figure 
one can identify at least four different fans associated with the transitions between a = 0, 1 
hole and /3 = 0, 1 electron states. However, in the region of strong static fields considered here, 
the majority of these transitions are weak and the whole spectrum is dominated by the vertical 
L = transition between the ground hole and electron states. Note a complicated structure 



''The superlattice parameters are Vb = 0.0632 eV (Vb — —0.0368 eV) for the electron (hole) potential barrier, and 
m* = 0.067me (m* — 0.45me) for effective electron (hole) mass. These parameters correspond to the value of the 
scaled "electron" and "hole" Planck constants h = 3.28 and h = 1.64, respectively. 
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Figure 4.7. Grey-scaled map of the one- dimensional absorption spectra ( 4-21 ) a function of 
the static field F and photon energy hu. 



of the main line resembling a broken feather. Recalling the results of Sec. (see Fig. 
this structure originates from avoided crossings between the (ground) under-barrier and (first) 
above-barrier electron resonances. Such a "broken feather" structure was well observed in the 
cited experiment |121]. 



Chapter 5 

Quasienergy Wannier-Stark states 



In the following chapters we investigate Wannier-Stark ladders in combined ac and dc fields. 
Then the Hamiltonian of the system is 

(5.1) 



H = Y + + Fx + F^x cos{ujt) , 



or, as described in Sec. equivalently given by 

2 

H = ^ + V[x + ecos{ujt)]+ Fx, e = F^/u^ . (5.2) 

Depending on the particular analytical approach we shall use either of these two forms. Let us 
also note that the Hamiltonian (|5.2| ) can be generalized to include the case of arbitrary space- 
and time-periodic potential V{x, t) = V{x + 27r, t) = V{x, t + T^). 

5.1. Single-band quasienergy spectrum 

For time-dependent potentials the period of the potential sets an additional time scale. In 
order to define a Floquet-Bloch operator with properties similar to the time-independent case, we 
have the restriction that the period of the potential and the Bloch time Tg are commensurate, 
i.e. 

pT^ = qTB = T . (5.3) 
In this case the Floquet operator U (T) over the common period T can be presented as 



U{T) 



-iqx 



U{T) , U{T) = e5q3 



dt 



{p - Ftf 



V{x,t) 



(5.4) 



(compare with Eqs. ( |2.10D - p.llD ). Consequently the eigenstates of U{T), 

U{T) $(x; A, k) = \^{x; A, k) , A = exp(-i£;r/;i) , (5.5) 

can be chosen to be the Bloch-like states ||177| , |178|| , i.e. <I>(x + 27r; A, k) = e^^'^'^$(a;; A, k). Due 
to the time-periodicity of the potential, V{x,t + T^) = V{x,t), we have the relation 

U{T) = U{T^Y = fexp(-ixg/ij) ^(r^)l'' . (5.6) 
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As a direct consequence of this relation, the states ^{x; A, k) with the quasimomentum k — r/p 
(r = 0, 1, ... ,p — 1) are Floquet states with the same quasienergy. In terms of the operator 
U^'^\T) = exp(— iKx)?7(r) exp(iKx) this means that the operators U^'^\T) are unitarily equiv- 
alent for these values of the quasimomentum]^ Therefore, the Brillouin zone of the Floquet 
operator U(T) is p-fold degenerate. In the next section we introduce the resonance Wannier- 
Bloch functions ^a,K{x) which satisfy the eigenvalue equation (|5.5| ) with the Siegert (i.e. purely 
outgoing wave) boundary condition and correspond to the complex energy ^^(k). Then the p- 
fold degeneracy of the Brillouin zone just means that the dispersion relation iSq(k) is a periodic 
function of the quasimomentum with period given by p. 

It should be noted that the Wannier-Bloch functions ^{x; A, k) (hermitian boundary condition) 
or $a,K(x) (Siegert boundary condition) are not the quasienergy functions of the system because 
the latter, by definition, are the eigenfunctions of the evolution operator U{T^) over the period 
of the driving force. However, the quasienergy functions can be expressed in terms of the 
Wannier-Bloch functions as 



^ ^ '" 27rn 

-i r 

p 



exp 

^ r=0 



'^a,K+r/pix) ■ (5.7) 



Equation ( ^.7] ) is the discrete analogue of the relation (|2.40| ) between the Wannier-Bloch and 
Wannier-Stark states in the case of pure dc field. Since the evolution operator U(T^) commutes 
with the translational operator over p lattice periods, the quasienergy states ^'^"^(a;) are the 
eigenfunctions of this shift operator. In particular, as easily deduced from Eq. (^.Tj), in the limit 
e ^ the function '^^J^{x) is a linear combination of every p-th state of the Wannier-Stark 
ladder (and altogether there are p different subladders). Thus, as well as the Wannier-Bloch 
states ^a,Kix), the eigenstates of U{T^) are extended states. Note that the Brillouin zone is 
reduced now by a factor p, i.e the quasimomentum is restricted to —l/2p < k < l/2p. On 
the other hand, as = T/p, the energy Brillouin zone is enlarged by this factor, i.e. the 
quasienergies take values in the interval < Re£^ < huj. Thus, if Sai^) is the complex band of 
the Floquet operator (|5.4|) , the complex quasienergies corresponding to the quasienergy states 
( IsTtD are 

4")(k) =f„(K) + , nio = 2TTF^. (5.8) 
p q 

In the remainder of this section we discuss the dispersion relation Sai^) for the quasienergy 
bands on the basis of the single-band model. It is understood, however, that the single-band 
approach can describe at its best only the real part E = HeS of the spectrum. 



In the single band analysis |54|, it is convenient to work in the representation (p. 11). Assuming 



that the two timescales are commensurate, the Houston functions ( 1.13| ) can be generalized to 



the Wannier-Bloch functions, which yields the following result for the quasienergy spectrum 

1 Ft 

Fain) = — I ea{n{t))At , K{t) = k — - sin(u;t) , (5.9) 

In this equation, as before, ea{K) is the Bloch spectrum of the field- free Hamiltonian Hq = 
p'^/2 + V{x) and K{t) is the solution of the classical equation of motion for the quasimomentum 

^We recall that in the case of pure dc field the operators U^'^\Tb) are unitarily equivalent for arbitrary k. 
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Figure 5.1. The band structure of the quasienergy spectrum originating from the ground (a = 0) 
Bloch hand as predicted by the single-band model. The parameters of the non- driven system are 
h = 3, F = 0.08 and the driving amplitude is e = 1. Only the rational values lo/ub = p/q with 
q < 7 are considered. The straight lines restricts the interval \E\ < hu corresponding to two 
(quasi)energy Brillouin zones. 



with initial value k. Expanding the Bloch dispersion relation into the Fourier series 

oo 

= ^^ea(i^) cos(27ri^K) (5.10) 

i/=0 

we obtain after some transformations 

oo ^qF 

//=0 

Thus, the dispersion relation for the quasienergies is given by the original Bloch dispersion 
relation with rescaled Fourier coefficients. For the low- lying bands, the coefficients ?a(i^) rapidly 
decrease with and for practical purpose it is enough to keep only two first terms in the sum 
over ^. 

Because the absolute value of the Bessel function is smaller than unity, the width of the 
quasienergy band is always smaller than the width of the parent Bloch band. In particular, 
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assuming €a{i^) ~ ea + (AQ,/2) cos(27rK) (as in the tight-binding approximation) and the simplest 
case of the resonant driving to = ujb (p = g = l), we have 

Ea{K) « e« + Ji (^) ^ cos(27r/t) . (5.12) 

As follows from this equation, the width of the quasienergy band approaches zero at zeros of 
the Bessel function Ji{z). This phenomenon is often referred to in the literature as a dynamical 
band suppression in combined ac-dc fields |4^5^ ^ A similar behavior in the case of a pure ac 



field was predicted in |49,^| and experimentally observed in |13J]. 

Let us finally discuss the case of an irrational ratio of the Bloch and the driving frequency, 
7 = lu/lob- We can successively approximate the irrational 7 by rational numbers Pj/qj, which 
are the j-th approximants of a continued fraction expansion of 7. Then, as for a typical 7 
both Pj,qj 00, the bandwidth of this approximation exponentially decreases to zero and the 



quasienergy spectrum turns into a discrete point spectrum ||5^. This is illustrated by Fig. |5J. 
where the band structure of the quasienergy spectrum (5^), calculated on the basis of Eq. ( 5.11 ), 
is presented for a = and constant value of driving amplitude e = F^lu'^. (The parameters of 
the non-driven system with V{x) = cosx are h = 3 and F = 0.08.) Note that the quasienergy 
bands have a noticeable width only for integer values of p. 

It is an appropriate place here to note the similarity between the quasienergy spectrum of a 
driven Wannier-Stark system and the energy spectrum of a Bloch electron in a constant magnetic 
field. The latter is known to depend on the so-called magnetic matching ratio 

where d is the lattice period. The spectrum of the ground state energies as a function of /? forms 
the famous Hofstadter butterfiy |18C]. In particular, for rational control parameter /? = p/q 



the number of distinct energy bands in the spectrum is given by the denominator q. Note that 
the magnetic matching ratio can be interpreted as ratio of two timescales, one of which is the 
time (i^m/2'Kh a particle with momentum 2iTh/d needs to cross the fundamental period d, and 
the other is the period eB/mc of the cyclotron motion.^ Similar, the driven Wannier-Stark 
system has two intrinsic timescales and the structure of the quasienergy spectrum depends on 
the control parameter 7 = Tb/T^ = fiio/edF, which is often referred to as the electric matching 
ratio. 



5.2. S-matrix for time-dependent potentials 



Provided the condition (5.3) is satisfied, the definition of a scattering matrix closely follows 
that of Sec. 2.2. Thus we begin with the matrix form of the eigenvalue equation ( |5.5| ), which 
reads 



(5.14) 



^Actually this phenomenon (although under a different name) was known earlier [179 
^This remark is ascribed to F. Bloch. 
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(To simplify the formulas we shall omit the quasimomentum index in what follows.) Comparing 



this equation with Eq. (2.18|), we note that index of the matrix U is now shifted by q. Because 



of this, we have q different asymptotic solutions, which should be matched to each other. Using 
the terminology of the common scattering theory we shall call these solution the channels. 

It is worth to stress the difference in the notion of decay channels introduced above and the 
notion of decay channels in the problem of above threshold ionization (a quantum particle in a 
single potential well subject to a time-periodic perturbation) |181| ] . In the latter case there is a 
well defined zero energy in the problem (e.g., a ground state of the system). Then the periodic 
driving originates a ladder of quasienergy resonances separated by quanta hco of the external 
field and, thus, the number of the corresponding decay channels is infinite. In the Wannier-Stark 
system, however, the ladder induced by the periodic driving (let us first discuss the simplest 
case p = q = 1) coincides with the original Wannier-Stark ladder. In this sense the driving does 
not introduce new decay channels. These new channels appear only when the induced ladder 
does not coincide with the original ladder. Moreover, in the commensurate case lo/lvb = p/q 
(because of the partial coincidence of the ladders) their number remains finite. With this remark 
reserved we proceed further. 

As before, we decompose the vector Gs into three parts, i.e. contains all coefficients with 
n > N and ^ all coefficients with n < —N — q. The third part, G^p, contains all remaining 
coefficients with — — q < n < N . The coefficients of and ^ are defined recursively, 

Gsim) = {X/um)Gsim- q) for m > A , (5.15) 

Gs{m- q) = {um/X)Gs{m) form<-A, (5.16) 

where Um = exp(i/i^[(K -\- m — q)^ — {k + ■mf]/QF). Let W be the matrix U truncated to the 
size (2 A -|- 1) X (2 A -|- 1), and, furthermore, let Om,n be an m x n matrix of zeros. With the 
help of the definition 




02N+l,q 

the equation for G^g' reads 



(5.17) 



[ UN+qGs{N + q) \ 

(5.18) 



UN+iGs{N + l) 



{Bn - Al) Cf = - 

\ 02N+1,1 J 

The right hand side of the last equation contains q subsequent terms Gs{m) and therefore 
contributions from the q different incoming asymptotes. However, we can treat the different 
incoming channels separately, because the sum of solutions for different inhomogeneities yields 
a solution of the equation with the summed inhomogeneity. Thus, let us rewrite ( 5.1^ ) in a way 
that separates the incoming channels. We define the matrices e'^ and Cq as 

e''=(^^'''' ), e, = (Og,27V+i,l,,,) , (5.19) 

\ 02N+l,q J 



49 



where Iq denotes a unit matrix of size qx q. Furthermore, we define the matrix u'^ as a diagonal 
q X q matrix u'' with the diagonal 

diag('u'') = (uN+q, . . . , UN+i) (5.20) 

and finally the column vectors C and Gq with the entries G{N + q), . . . , G{N + 1) and G{—N — 
1), . . . , G{—N—q), respectively. With the help of these definitions the right hand side of equation 
( |5.18| ) reads e'^u'^G'^, which directly leads to the following relation between the coefficients of the 
incoming and the outgoing channels 

Gq = eq[BN - Xt]-\iu'^G'^ . (5.21) 

In the S-matrix formula we additionally need to include the influence of the free states, which 
are again discrete versions of Airy functions. Thus, with the help of two additional diagonal 
matrices, a'^{E,N) and aq{E,N), which contain the contributions of the free solutions, 

diag(a'?) = {Go{N + ?),..., Go{N + 1)) , diag(ag) = {Go{-N - 1), ... , Go{-N - q)) (5.22) 

with Go{m) = exp (^ih?[K + m]^ /6F — iE[K + m]/F^, we define the q x q S-matrix 

S{E)= lim aq^ eq[BN - Xir^ e'^u'^ a'^ . (5.23) 

N^oo 

It can be proved that the matrix ( ^.23| ) is unitary by construction, i.e. {E)S{E) = 1. 
Based on Eq. ( [5.23| ), the equation for the resonance wave functions has the form 

{Bn-X1)gP=0. (5.24) 

In fact, as follows from the explicit form of the matrix B^, the first q elements of the eigenvector 
are zero and, according to Eq. ( 5.16| ), G^^^ = 0. Thus, the solution of Eq. ( 5.24 ) satisfies the 
resonance-like boundary condition of empty incoming channels. The corresponding energies are 
given hy £ = ihlnX/T and actually depend on k, which enters all equations displayed above as 
a parameter. 

To conclude this section, we generalize the equation for the Wigner delay time. The general- 
ization of ( |2.30| ) for systems with q decay channels reads 

ih dln[detSiE)] 

J dE ' ^^-^^^ 
or, equivalently 

r = ^Tr(f), ?=-ihSHE)^^. (5.26) 



where r is the so-called Smith matrix |182|. Along with the Wigner delay time, in the random 
matrix theory of chaotic scattering (see chapter 7) the notion of partial delay times, which are 
the eigenvalues of the Smith matrix, and one-channel delay times, which are the the diagonal 
elements of the Smith matrix, appear. The sum of the partial or one-channel delay times 
obviously yields the Wigner delay time. 
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Figure 5.2. T/ie rea/ (left panel) and the imaginary (right panel) parts of the ground quasienergy 
resonances as function of the driving frequency cu. The dashed line interpolates the average 
values Eq{uj) and f o(w) obtained for the rational values u!/ujb = p/<l with q < 7 and p < 98. 
The "error bars" mark the bandwidths Aq'^{u!) and AQ™(a;). The system parameters are h = 2, 
F = 0.061 and e = 0.08. 



5.3. Complex quasienergy spectrum 

Using the scattering matrix approach of the preceding section we can calculate the complex 
quasienergy spectrum of the Wannier-Stark system for arbitrary values of the parameters. In 
this chapter, however, we confine ourselves to the perturbation regime of small e and relatively 
large values of the scaled Plank constant h. The opposite case of large e and small h will be 
considered in chapter 7. 

We begin with the analysis of the real part of the spectrum, E = Ref . Recalling the results 



of Sec. |5jJ the real part of the quasienergy spectrum is expected to obey 

X 2TTFn A?-° 
EI^>{k) =E^ + + cos{2ttpk) , n = 0, . . . ,p ■ 



1 



(5.27) 



The left panel in Fig. ^ shows the mean position of the ground quasienergy bands (dots) and the 
band widths (marked as error bars) calculated for some rational values of the driving frequency 
u! (only the bands with n = are shown). The parameters of the non-driven system with 
V{x) = cosx are h = 2 and F = 0.061. For these parameters the widths of two first resonances 
are Tq = 1.24 • 10~^ and Fi = 1.30 • 10~^. The distance between the real parts of the resonances 
is El — Eq = 3.78AHujb- It is seen in the figure that the band widths Aq'^ = AQ''^(a;) are large 
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only for lo = pcoB, in qualitative agreement with the estimate ( 5.11 ). We would like to stress, 
however, that estimate ( [5.11 ) is obtained within the single-band approximation and, because 
of this, the actual bandwidths deviate from this dependence. (We shall discuss the conditions 
of validity of Eq. (5.11) later on in Sec. ^^.) The second deviation from the predictions of 
single-band model is the dependence of the mean quasienergy band position Eq on lo. As shown 
below, this dependence reflects the presence of the other quasienergy states, originating from 
the higher (a > 0) Bloch bands. Let us also note that the mean position Ea = Ea{uj) is, unlike 
the band width A^^ = A^*'(a;), a continuous function of the frequency. 

The right panel in Fig. |5.2| shows the imaginary part F = — 2ImiS of the quasienergy spectrum. 
In the perturbation regime e 



a behavior similar to ( 5.27 ) 



F„(k) PS F„ 



A 



Im 



cos(27rK) 



(5.28) 



is observed. It should be noted that the smooth function Fq = Fo(a') approximating the mean 



values of the bands is nothing else as the induced decay rate discussed in Sec. 4.1. In fact, an 



arbitrary initial state of the system (which was assumed to be the ground Wannier-Stark state 
^o,/(2^) ill Sec. |4.1| ) can be expanded in the basis of the quasienergy states ^q"'*(x) as 

(5.29) 



^{t = 0) = ^c„,„*W , c„(n) = (M/(r)|*(t = 0)) . 



(Here we assume that uj/ujb is an irrational number and, therefore, the quasienergy functions 
are localized function with discrete spectrum.) During the time evolution the coefficients Ca,nit) 
decay as exp{—Tat/2h). Since F^ > Fq (a > 0), the projection of the wave function back to the 
initial state decays (after a short transient) exponentially with an increment given by ro(u;)/2. 
This is the underlying argument of our numerical method of calculating the decay spectrum of 
the system. Namely, to obtain the decay spectrum discussed in chapter 4 we calculated the mean 
imaginary values of the quasienergy bands for a number of rational uj/ub and then interpolate 
them for an arbitrary u). 

Let us now discuss the w-dependence of the smooth functions Eq{oj), Tq{uj). Because we ana- 
lyze the case of weak driving, these functions can be obtained by using perturbation theory. In 
fact, assuming again an irrational value of w/cjb, the zero order approximation of the most sta- 
ble quasienergy function is the ground Wannier-Stark state ^'o,n(a^)- According to the common 
perturbation theory, the first order correction is 



(n) 



0,n 



(^0,n|2;|^ 



a,l) 



il ± 



Correspondingly, the second oder correction to the energy is 



(n) 



£o,n + 



EE 



n 



£a — £o + {I — n)hjjJB ± ^ 



(5.30) 



(5.31) 



In Eq. (5.31) we used the notation (4.3) for the squared dipole matrix elements and took into 
account that the energies of the Wannier-Stark states form the ladder E^^i = £a+l^B- Equation 
( |5.31| ) is illustrated in Fig. 5.3, where the real (left panel) and imaginary (right panel) parts of the 
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Figure 5.3. Corrections to the ground state energy from Fig. \5.l\ (dashed line) compared to 
approximations based on equation ( [5. 3^ ) (solid line). The left panel shows the real part of the 
ground state energy, the right panel the imaginary part. 



quasienergy calculated on the basis of this equation (solid line) are compared to the numerical 
data of Fig. (dots, interpolated by a dashed line). For small (relative to ujb) frequencies 
both curves coincide almost perfectly, but deviate for large oj. This deviation can be attributed 
to the slow convergence of the perturbation series over a in the high-frequency region. (For the 
presented results, the upper limit for the sum over the Wannier-Stark ladders is taken as a = 3.) 



The concluding remark of this section concerns the relation between Eq. (O) [i.e. the imag- 
inary part of Eq. (|5.31| )] and the "e-version" of Eq. ( |4.2| ) used to analyze the decay spectrum of 
atoms in optical lattices in Sec. O. The difference is the use of the squared matrix elements 
( 4.171 ) instead of the squared dipole matrix elements ( |4.3| ). However, recalling the relation 
e = F^juP' and relation ( [4.19 ), this difference can be shown to be within the accuracy of the 
second order perturbation theory. The advantage of the e-version over the i<lj-version is a better 
convergence in the high-frequency region. 



5.4. Perturbation theory for rational frequencies 

Discussing the perturbation approach in the previous section we excluded the case of rational 
ratio of the driving and Bloch frequencies. Let us now turn to it. To be concrete, we restrict 
ourselves by the simplest but important case lo = oje- In this case the periodic driving couples the 
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Wannier-Stark states belonging to the same Wannier-Stark ladder and, therefore, the extended 
Wannier-Bloch function ^a,K is an appropriate zero order approximation to the quasienergy 
function. 

As described in the beginning of this chapter, the complex quasienergies of the system are 
found by solving the eigenvalue equation 

U{T^)^a,^ = exp[-i£a{K)T^/h] . (5.32) 

Let us approximate the Hamiltonian (^]^) by the first order of the Taylor expansion in e, H ^ 
Hw — e sin(x) cos(wt) (here H]y is the Wannier-Stark Hamiltonian ( |2.2| ) and V{x) = cosx is 
assumed for simplicity). Then we can calculate the effect of the periodic driving in the interaction 
representation of the Schrodinger equation. Explicitly, we get 

U{T^)^Ue{T^)Uw{T^) (5.33) 
where the operator Uir(T^) reads 

Ue{T^) = eip ^ " dt cos(L^t) U^it) sin(x) Uw{t)^ , (5.34) 

and the operator Uwii) is the evolution operator for the unperturbed system. According to 
common perturbation theory, the first order correction is given by the diagonal elements of the 
operator U£{T^), 

exp[-iA£^{K)T^/h] = Ue{T^) . (5.35) 

Let us approximate this formula further. Expanding the operator exponent in a series in e and 
keeping only the first term, the correction to the quasienergy reads 

AEa{K) = -7^ dt cos(wt) Ul^^it) sin{x) Uwit) \^a,K,) ■ (5.36) 

Using the solution Uw(t)^a,K = exp(— i<Sat//i) ^a,K~Ft/h substituting dt/T^ = —dn, 
Eq. ( ^.36| ) takes the form 

A£a{K) = -e dK cos{2ttk') {^a,K+K'\sin{x)\^a,^z+^z') ■ (5.37) 
Jo 

Finally, using the symmetry property of the Wannier-Bloch function, the integral ( 5.37| ) can be 
presented in the form 

A£aiK) = " ^ " cos(27rK) , (5.38) 



where 

A^*^ + iAjj^ = -2e /" dKcos(27rK) sin(x) I (5.39) 

Jo 

(The special notation for the band width stresses that the integral on the right hand side of 
Eq. ( p.39| ) is a complex number.) Thus, a weak periodic driving removes the degeneracy of the 
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Wannier-Bloch bands which then gain a finite width. Moreover, there are corrections both to 
the real and imaginary part of the quasienergy. 

In conclusion, let us briefly discuss the relation between the formulas ( 5.3S| ), (15.39 ) and the 
tight-binding result ( 5.12| ). As was stated many times, the single-band model neglects the 
interband tunneling, which is justified in the limit F — > 0. In this limit the quasienergy band 
width can be estimated as 

47ri^ZX„ . ^^^^^^ 



A 



Re 



7? 

A — 



^2 ~" F 

where Aq is the width of the Bloch band. Indeed, using Eq. dilsD , the band width in Eq. ( |5^ ) 
can be expressed in terms of the dipole matrix elements as 



A^^ + iA: 



Im 



{2TiFf 



(5.41) 



Then, using the tight-binding approximation ( |1 . 1 1 ) for the resonance Wannier-Stark states 
^a,i{x)^ we obtain the estimate ( 5.401 ). (Alternatively, we can approximate Xa,K{x) in Eq. ( 4.14| ) 
by the periodic part of the Bloch function.) It is seen, that the estimate ( 5.40 ) coincides with 
Eq. ( ^.12 ) in the limit F^^/F 0. We would like to stress, however, that the actual perturbation 
parameter of the problem is e ~ F^jF"^ and not F^jF ^ as it could be naively expected on the 
basis of the tight-binding model. 



5.5. Selective decay 



This section serves as an illustration to the perturbation theory of Sec. 5.4 and discusses some 
important limitations of the perturbation approach. In order to reduce the number of relevant 
resonance states, we choose the parameters of the unperturbed system as = 2, F = 0.08. In this 
case we have to take into account mainly two resonances with energies = 9.42-10~2— i5.60T0~'^ 
and £i = 4.18 • 10"^ — i8.81 • lO"^. All other resonances are very unstable and approximately 
do not infiuence the results. The frequency of the time-periodic perturbation is given by w = 
27rF/^« 0.251. 



Figure 5^ shows the real parts of the quasienergies of the two most stable Wannier-Bloch 
resonances for different amplitudes e. In panel (a) we have the unperturbed case with fiat bands. 
When the ac driving is added, the dispersion relation of the ground band is well described by the 
theoretical cosine dependence. The first excited band follows this relation only up to e = 0.2. 
If the amplitude is increased further, deviations from the cosine appear, and for e = 1 other 
effects strongly infiuence the band (note that in this case F^^/F 0.79, thus we are still far away 
from the parameter range where the tight-binding model predicts dynamical band suppression) . 
Furthermore, for e > 0.2 the bands cross, and then we cannot neglect their interaction. 

Next we investigate the bandwidth, i.e. the difference between the extrema of the real parts 
of the quasienergies, AE = E{k = 1/2) — E{k = 0). Figure 5.5 shows the width of the two 
most stable bands as a function of the amplitude e for three different field strengths F = 0.02, 
F = 0.04 and F = 0.08. It is seen that in all cases the bandwidth grows approximately linearly. 
Again, the agreement is much better for the ground band; for the first excited band one observes 
an oscillation around the linear growth. Note that the slope is proportional to F as expected 
on the basis of the perturbation theory [see Eq. ( 5.4C| )]. 
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Figure 5.4. Real part of the two most stable (quasi-) energy Wannier-Bloch bands for e = (a), 
£ = 0.2 (b), e = 0.4 (c) and e = 1 (d). The other system parameters are h = 2, F = 0.08 and 
oj = 27rF/;i« 0.251. 




Figure 5.5. Band width of the two most stable Wannier-Bloch bands as a function of the per- 
turbation parameter e for F = 0.08 (solid line), F = 0.04 (dashed line) and F = 0.02 (dotted 
line). 
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Figure 5.6. Width (decay rate) of the ground state as a function of the Bloch index k for the 
cases studied in Fig. \5.4 - The width is normalized with respect to the width at e = 0. 




Figure 5.7. Decay rate of the most stable Wannier-Bloch state at k = (lower family of curves) 
and K = ±1/2 (upper family of curves) for the same parameters as in Fig. |5.4 i-e. F = 0.08 
(solid line), F = 0.04 (dashed line) and F = 0.02 (dotted line). The rate is normalized against 
the decay rate at e = 0. 
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We proceed with the analysis of the imaginary part of the quasienergy spectrum. Figure 5.6 



shows the width of the ground state as a function of the Bloch index for the parameters of 



Fig. |5.4 For e = the band is flat as predicted from the theory. For e = 0.2 the width can be 
approximated by a cosine, however, the mean is shifted to approximately twice the unperturbed 
width. If we further increase e, additional structures appear. In comparison with Fig. ^.4| , we 
see that the bandwidth is increased where the (real part of the quasienergy of the) ground band 
crosses the first excited band. Therefore, we can clearly assign the increase of the width to the 
band crossings. Recalling the results of chapter 3, we again observe effects of resonant tunneling, 
now as a function of the quasimomentum. As shown in Ref. [156,164], the two-state model of 



Sec. can be adopted to the present case and yields good correspondence to the numerical 
data. 



In Fig. p.6| we can see that the perturbation can both increase and decrease the width and 
thus the rate of decay of the quasienergy states. In the case considered, for small e the decay is 
enhanced at the edges of the Brillouin zone and suppressed in its center Let us therefore take 
these two quasimomenta to further investigate the dependence on the perturbation parameter e. 
The results of a calculation of the widths at k = and fi; = ibl/2asa function of the amplitude 
e are shown in figure [S.T] . For small e the dependence is nearly linear, but for larger values it is 
highly nontrivial. In particular, we would like to draw the attention to the behavior of the solid 
line at e ~ 0.4 and of the dashed line at e ~ 0.54. Here the decay rate is suppressed by more 
than a factor 10^! This tremendous decrease of the decay rate has enormous consequences on 
the global dynamics. For example, let us initially take the most stable Wannier-Stark state and 
then add the ac driving. Then the survival probability is given by 

P(t) = y"'^^'dKexp(-^^) . (5.42) 

If we approach the critical value Ecr, the decay is suppressed and asymptotically 

Pit) - exp(-r^i„ t/h) , (5.43) 

where Fmin is the minimal decay rate. Let us also note another property. Since the decay rate of 
the quasienergy states depend on the quasimomentum, after some time only the contributions 
with quasimomentum around the value with the smallest decay rate will survive. In what follows 
we shall refer to this phenomenon as the selective decay of the quasienergy states. Some physical 
consequence of this phenomenon are discussed in the next chapter. 



*The regions of enhanced and suppressed decay depend on the difference between the phase of the driving force 
and the phase of the Bloch oscillation. For example, the change of cos{ujt) in the Hamiltonian ( ^.2| ) to sm{Ljt) 
shifts the displayed dispersion relation by a quarter of the Brillouin zone. 



Chapter 6 

Wave packet dynamics 



In this section we address the question of the time evolution of an initially localized wave packet. 
Usuahy this problem is analyzed by simulating the wave packet dynamics on the basis of the 
time-dependent Schrodinger equation. However, this numerical approach is very time consuming 
and has an upper limit for the times considered. In what follows we describe the evolution of the 
wave packet in terms of the resonance states. Besides tremendous decrease of the computational 
efforts, the latter approach also gives additional insight into the decay process of the Wannier- 
Stark states. 



6.1. Expansion over resonance states 

A direct expansion of a localized state in terms of resonances yields inappropriate results 
because in the negative x-direction the resonance states extend to infinity. Therefore the de- 
scription needs to be modified to take into account the finite extension of the initial state. 
Recently this problem was analyzed for decaying quantum systems with a finite range potential 



Let us adopt the approach of |184| to describe the evolution of the wave packet in momen- 
tum space. In this approach, the wave function ^{k, t) is expressed in terms of the stationary 
scattering states ^s{k',E): 



/oo / 17^ \ 

dE f{E) ^s{k-, E) exp f -i — j , (6.1) 

where f{E) = {^l s{k]E)\tjj{k,Q)). [We recall that the states ^s{k]E) are normalized to 5- 
function: {"^ s{k\ E')\^ s{k\ E)) = 6{E — E').] We are mainly interested in the properties of 
the decay tail at A; — > — cxd. In this region the scattering states can be approximated by their 
asymptotic form [see Eq. 



^2^.3 Fk\ , N 

hm *±(A;; E) = g±{E) exp i -- - i — , g±{E) = e±''^(^) . (6.2) 



^However, this problem was already addressed in textbooks as, e.g., [185| 
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Substituting this asymptotic form into Eq. (^]^) we have 

tl;{k,t) = ex.p(i^^\ G-(k + ^\ , k <€. , (6.3) 



. 6F / V h 

where 



oo 



G.ik)=l dEl^exp(-i^) . (6.4) 



oo 



If the initial wave function ilj{k, 0) has a finite support, f{E) is an entire function in the complex 
plane. Then the function f{E)/g^{E) has simple poles at zeros of g+{E), i.e at the poles of 
the scattering matrix S{E) = g-{E)/g+{E). This property suggests to evaluate the integral 
( |6.4D with the help of the residuum theorem. Without knowing the explicit form of the function 
f{E)/g+{E) we have to make some assumptions on its asymptotic behavior in order to proceed 
further. In particular, if we assume that the function f{E)/g^{E) does not influence the be- 
havior of the integrand at infinity, the integral yields a sum over the residua located within the 
appropriate contour. Explicitly, for k > the contour should be closed in the lower half of the 
complex plane, for A: < it contains the upper half. Since all poles of the scattering matrix are 
located in the lower half of the complex plane, we get 

G_ (fc) = 2m e (k) Y^b,exp(-i^Y (6.5) 

where Q{k) is the Heaviside function, the bi, are the residua of f{E)/g^{E) at the poles, and 
v = {a, I}. Inserting this result in ( |6.3D yields 

u c +\i 

(6.6) 



V'(fc, t) = e {hk + Ft) ^ c. exp [i - ^ - ^) 



with Cu = 2TTibu. The terms of the sum are actually proportional to the asymptotic form of the 
resonance wave functions ^u{k,t). Thus, we can equivalently present the wave function as 



Qihk + Ft)Y,Cue^p(^-ijpj^u{k) . (6.7) 



Therefore, in the Stark case we can describe the evolution of an initial state by a superposi- 
tion of resonances, where we take into account the space-time decay process in the prefactor 
Q{hk + Ft). This factor truncates the wave function at the momentum hk = —Ft, i.e. only 
momenta with hk > —Ft contribute. With increasing time, the wave function extends to smaller 
momenta, where the edge moves according to the classical equation of motion. 

It should be noted that the location of the edge reflects the assumption on the behavior 
at infinity we made in order to explicitly evaluate the integral. For example, the function 



f{E)/g^{E) may contain an additional exponential factor exp(ia£') (see the example in [184|). 
Though this factor does not influence the poles, it nevertheless influences the argument of the 
Heaviside function. In fact, in a realistic situation the edge will be shifted, because the truncation 
edge at i = has to reflect the extension of the initial state in momentum space. We take this 
into account by replacing hk in the argument of the Heaviside function by h{k + ko), where 
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ko describes the extension of the initial state in the negative fc-direction. Furthermore, if the 
initial state does not have a compact support but a tail in the negative momentum direction, 
the edge will be smoothed and deformed. However, the qualitative behavior remains unchanged: 
the prefactor is approximately constant for positive arguments of the Heaviside function, and 
it approximately vanishes for negative arguments. Therefore, we take the Heaviside description 
as a reasonable approximation to the real situation. Let us also note that the wave function 
constructed in this way can be normalized. Indeed, in the positive momentum direction the 
resonances decrease stronger than exponentially, and in the negative direction the wave function 
is truncated. 

Now we discuss the dynamics of the wave packet in coordinate space. If we are interested 
in the asymptotic behavior for x <^ 0, the wave function t/j^x, t) can be found by a Fourier 



transform of the asymptotic form of Eq. (6.7): 



/oo 
dke[h{k + ko) + Ft] ^ Cu exp 
-oo 



2l3 



K'k 



Eyk 



h 



+ kx 



(6.8) 



Let us evaluate the integral in the stationary phase approximation. The equation for the sta- 
tionary phases reads 



dA? 



6F 



k£u 



+ kx 







(6.9) 



Neglecting the imaginary part of the energy Si, ^, the stationary phase condition is just the energy 
conservation, and the stationary points are the classical momenta hk^, = ^J2{E^ — Fx) = py{x). 
If Pv{x) <C —hk{) — Ft, the prefactor is zero and the integral vanishes. On the other hand, if 
Pu{x) ^ —hko — Ft, the integral of the contribution of the ly-th resonance yields approximately 



exp 




exp 



.pI{x) T^p^ix) 



3hF 



2hF 



(6.10) 



which is just the asymptotic form of the Wannier-Stark state in the coordinate representation. 
The critical point is Pu{x) = —hkQ — Ft, where the approximation breaks down because the 
Heaviside function is not a slowly varying function at this point. Actually, in the vicinity of this 
point the integral interpolates between the other two possibilities. Let us skip a more detailed 
analysis here and roughly describe the transition between both regimes by a Heaviside function 
of the argument Pi,{x) + hko + Ft, or, equivalently, of the argument x + F{t + to)^/2 — E^/F, 
where to = hk^/F. Then, replacing the contribution ( |6.10 ) by '^u{x,t), we get 



i^{x,t) = 9 



X + 



F{t + tof 



F 



exp 



h 



^^uix). 



(6.11) 



In comparison to equation (|6.7] ) there are two differences. First, in coordinate space the trunca- 
tion depends on the energy of the resonances. Furthermore, the edges of the different contribu- 
tions move with a quadratic time dependence, which reflects the classical (accelerated) motion 
in a constant external field. 



^More precisely, we treat the exponential of the imaginary part as a slowly varying function. 
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6.2. Pulse output from Wannier-Stark systems 

Let us consider the dynamics of a coherent superposition of the Wannier-Stark resonances 
belonging to a particular Wannier-Stark ladder 

V'lfc,*) = X^Qexp ^-i^^^ ^'„^i(/c) , q = -^exp^-^^ . (6.12) 

(To shorten the notation, we skip here the truncation by the Heaviside function because the 
truncation does not influence the properties which we are going to discuss). This problem, as 
will be seen later on, is directly related to the experiment |126(| , where a coherent pulse output 
of cold atoms was observed. Based on this phenomenon, a possibility of constructing an atomic 
laser is currently discussed in the literature. 

According to Eq. ( p. 391) the Wannier-Stark states belonging to the same ladder are related by 
^a,i{k) = exp(—i2TTlk) ^a,oif^) cind Sa^i = £a + 27rlF. Combining this phase relation with the 
different phases due to the time evolution, the time evolution of the superposition is given by 

'^{k,t) = ^a,o{f^,t)C + , (7(A:) = ^ Q exp(-i27r/A;) . (6.13) 

where ^a,o{k,t) = exp(— i<SQ,t/^)^a,o(^)- Thus, the time evolution of the superposition is given 
by the time-evolved wave function at the mean energy, ^a^Q{k,t), times the discrete Fourier 
transform C{k) of the amplitudes q, which is taken at the momenta k + Ft/h. Since the 
function C{k) is periodic in momentum space the factor C{k + Ft/h) is also periodic in time 
with the period h/F = Tb- In what follows we shall refer to the function C{k) as amplitude 
modulation factor. In the considered case q ~ exp(— Z^/cr^) the amplitude modulation factor is 
obviously a periodic train of Gaussians with the width a~^. 

We turn to the coordinate representation. Following the derivation of the preceding section, 
the wave function ip{x, t) can be shown to obey 

i;{x, t) = ^^,oix, t)cUr + ^) , (6.14) 



where, as before, p{x) = ^j2[Ea — Fx) is the classical momentum. Because the function C{k) 
has peaks at integer values of the arguments, the function C^t/Ts +p{x)/h) has peaks at the 
coordinates 
p 

X = xq - — {t + mTB)"^ , (6.15) 

where xq = Ea/F is the classical turning point. Thus, as a function of time, the peaks accelerate 
according to the classical equation of motion of a free particle subject to a constant electric 
field. Additionally, the peaks broaden linearly with increasing time (or with increasing m). It is 



straightforward to combine the result ( 6.14 ) with the result of the previous sections. Generally, 
we have to truncate the wave front approximately at the coordinate x = —F{t + to)^/2. 

Figure |6.1| shows the evolution of the superposition of the ground Wannier-Stark resonances for 
V{x) = cosx. The system parameters are h = 3.3806, F = 0.0661 and a = 15, which correspond 
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Figure 6.1. Space-time decay of the wave packet for the parameters of the experiment \l2t ]. 
From top to bottom, the panels correspond to t = 3 ms, 5 ms, 7 ms and 10 ms, respectively. 



to the setting of the experiment [126|. The figure was calculated in the following way: First the 
ground / = Wannier-Stark state was calculated in the momentum representation. Then the 
wave function was multiplied with the amplitude modulation factor C{Ft/h + k) taken at the 
specified times and truncated according to equation (|6.7|). (We shifted the truncation edge by 
ko = 1/2 in order to avoid a truncation directly at the maxima. As mentioned in Sec. 6A this 



shift takes into account the finite extension of the initial state.) Finally the resulting function 
was Fourier transformed into coordinate space. The obtained result reproduces the findings of 



the experiment |126]. A series of pulses is formed which then accelerate according to the free 
motion. At a fixed value of the coordinate, the sequence is periodic in time (after the first pulse 
passed), up to an overall exponential decay which refiects the fact that every pulse takes away 
a certain amount of probability. 

A few words should be added about the validity of the one-particle approximation. In fact, in 
the experiment cited, the authors used a Bose-Einstein condensate of Rubidium atoms, uploaded 
in a vertically aligned optical lattice. Thus a description of the system with the help of the 
Gross-Pitaevskii equation^ 



j^ + Vo cos{2kLz) + Mgz + "^ntl^P 



(6.16) 



'A detailed introduction to the physics of Bose-Einstein condensates can be found in the review article ||l86l. 
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looks more appropriate. Equation ( 6.16 ) was studied numerically in Ref. [187-192]. It was 
found that for moderate densities of the condensate (realized in practice) the pulse formation 
is only slightly modified by the nonlinear term in the Gross-Pitaevskii equation. Thus the 
physics behind the experimentally observed phenomenon is provided by single-particle quantum 
mechanics and can be well understood in terms of Wannier-Stark resonance states. 

6.3. Atom laser mode-locking 

The crucial point for the existence of the pulse output in the Wannier-Stark system is the 



fixed phase relation between the probability amplitude q in Eq. ( 6.12 ). In the experiment 



[ |126| , this fixed ph ase relation was achieved by the self-interaction of the Bose condensate. In 



the following we show that one can prepare an appropriate initial state within single particle 
quantum mechanics. Explicitly, the statement is as follows. Take an arbitrary initial state (i.e. 
arbitrary ci ) and drive the system for a finite time Tint with the frequency matching the Bloch 
frequency ujb- If the driving amplitude is sufficiently large and the interaction time Tint is long 
enough, the initial state decays with a pulse output afterwards. 
The physics behind this effect is the selective decay of the quasienergy Wannier-Bloch state 



discussed in Sec. 5.5. Indeed, let the ^a,K{k) be the quasienergy states of the dc-ac Hamiltonian 



(5.2). Then we can expand the initial state tjj{k,Q) in this basis 

1/2 
/2 



/■1/2 

V'(A;,0) = V/ dAtc«(K)$«,«(A;), (6.17) 



where the Ca{n) are periodic functions of the quasimomentum. [In particular, assuming the 
adiabatic switching of the field, the initial condition ip{k., 0) = \I'o,/(^) will correspond to Cq,(k) = 
5a,o exp(iK/).] Note that the expansion ( |6.17| ) is also valid in the case of a pure dc field considered 
in the previous section, but we preferred there the alternative basis of the Wannier-Stark states 



[see Eq. (6.12)]. After N periods of driving the wave function reads 

^1/2 



^{k,NTB) 



l^^^dKCo(^) exp $^^^(^)^ (6.18) 



where we assumed that all quasienergy states, excluding the ground states a = 0, have decayed. 
Now the ac field is switched ofi^, and we take the final state ^{k^ NTs) as the initial state of the 
pure dc dynamics. Expanding it in the basis of the Wannier-Stark states yields ip{k,NTB) = 
co{k) exp{—i£o{k)NTB /h) ^ofi{k), where the functions co(k) and <So('^) are treated now as the 
periodic function of the momentum instead of the quasimomentum. Then 

ij{k, t>NTB) = ^.,o(A:, + , C{k) = co{k) exp f-i . (6.19) 



Comparing this result with Eq. ( |6.13 ), we notice that the prefactor co{k) exp{—i£Q{k)NTB /h) 



takes the role of the amplitude modulation factor C{k) of the new initial state. Let us discuss 
this factor in more detail. 



As shown in Sec. 5.4, for small e the dispersion relation of the complex quasienergy band is 



£o{i^) = "^o + (^0^*^/2 + iAo™/2) cos(27rK). Thus the absolute value of the amplitude modulation 
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Figure 6.2. Tail of the wave function after the system was driven for different times T = NTs- 
The parameters are h = 3.3806, F = 0.0661 and e = 0.1. 



factor is given by 
\C{k)\^ = |c(A:)pexp 



TqNTb A^^NTb 



h 



2h 



cos(27r/c) 



(6.20) 



If the interaction time Tint = NTs is large enough (and if c{k) is sufficiently smooth), the strong 
modulation of the exponential dominates the form of the amplitude modulation factor. Then 
the wave function is periodically peaked in momentum space. Of course, such a periodically 
peaked structure is also found for larger values of e where formula ( 5.3§| ) is no longer valid. In 
fact, due to the stronger modulation of ro(K), it appears even for short interaction times. 

The behavior of the wave function in coordinate space is additionally modified by the disper- 
sion due to the real parts of the quasienergies. If we approximate it by the cosine and again apply 
the stationary phase approximation in the Fourier transform of equation (|6.1S| ), the stationary 
points ks are solutions of the slightly modified equation ?i^A;^/2 + 7rAoReA^sin(27r/cs) = Eq — Fx. 
The implications are as follows. In coordinate space, the form of the peaks is changed compared 
to the dispersion-free case, in particular, the peaks can be broadened or narrowed. Note that 
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Figure 6.3. The same as in Fig. \6.^ but e = 1.5 



for small \ks\ there may be three instead of one stationary point for each branch of the square 
root. Then the wave function shows additional interferences due to the interaction of the three 
different contributions. However, for large |A;<j| (i.e. for x — > — oo), the dispersion only slightly 
influences the shape of the peaks. Thus, for large the shape of the peaks of the decay tail 
mainly reflects the function Tq{k), which provides a method to experimentally access of this 
function. 

To support the above analysis. Fig. |6.2| shows the decay tails which develop for a weak driving 
with £ = 0.1. In this numerical example we choose the ground Wannier-Stark resonance ^'0,0(3;) 
as the initial state (the other parameters are the same as in Fig. |6.lD and drive the system for 
different interaction time. (Explicitly, we calculated the resonances wave function \I'o,o(^)) mul- 
tiplied by the amplitude modulation factor exp(—i£Q{k)NTB/h), where the dispersion relation 
was calculated independently, and finally Fourier transformed to coordinate space.) After short 
interaction times, the tail is slightly modulated. For longer interaction times, the modulation 
depth increases and pulses develop, which finally are clearly separated. Note that, apart from 
effects due to the dispersion, we can decrease the width of the pulses by further increasing the 
interaction time, which provides a simple way to tune the width experimentally. 
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A crucial point of the weak driving regime is the long interaction time which is needed to 
generate well separated pulses. The relevant timescale is set by the most long-lived state from 
the ground band. For the case e = 0.1, the minimum width is rmin = 

7.214 • 10-3, which 

corresponds to a lifetime approximately 10 Tg. Thus, the interaction time is much longer than 
the lifetime of the most stable state. Consequently, a predominant part of the initial wave packet 
has already decayed before pulses are being formed. One can, however, surmount this problem 
by increasing the amplitude of the resonant driving. Figure 6.3 shows the decay tail for e = 1.5. 
Now the pulses develop after much shorter interaction times. For e = 1.5, the function Tq{k) 
has four minima, which are due to two crossings with higher excited Wannier-Stark ladders (see 
Sec. ^.4]) . Note that one can directly read off this property from the substructure of the pulses 
on the decay tail. In the lower panels of the figure one can also see the narrowing caused by 
the dispersion. In particular, the first peaks (counted from the right) strongly oscillate, which 
reflects the existence of three stationary points in this region. However, the last peaks have 
approximately the same shape, i.e. here the narrowing effect can be neglected. 



Chapter 7 

Chaotic scattering 



This chapter continues the analysis of Wannier-Stark system affected by an ac field. In chapters 



5-6 we have considered the case of a relatively large scaled Planck constant h (see Sec. 1.4) and 
relatively small values of driving amplitude e, where the perturbation approach can be applied to 
analyze the spectral and dynamical properties of the system. Now we turn to the region of small 
h ^ 0. In this region even a weak driving violates the condition of perturbation theory which 
roughly reads e/h < 1. On the other hand, a small h corresponds to the semiclassical region, 
where the classical mechanics can guide the quantum-mechanical analysis. It turns out (see next 
section) that the classical dynamics of the system ( |5.2D is typically chaotic. Then the question 
we address sounds as "What are the quantum manifestations of this chaotic dynamics?" . This 
question belongs to the list of problems considered by the modern branch of quantum mechanics 
known as Quantum Chaos (and actually can be considered as the definition of the field) ]194[|. 

A powerful tool of the theory of quantum chaos is the random matrix theory (RMT) |194- |197| . 
Its application is based on the conjecture that the spectral properties of a classically chaotic 
system are similar to those of a random matrix of the same (as the Hamiltonian) symmetry 
class. Recently a considerable progress has been made in nonhermitian random matrix theory, 
which aims at describing the properties of chaotic scattering systems [198-202|. In what follows 
we study the system (|5.2D from the point of view of nonhermitian random matrix theory. In par- 
ticular, we numerically calculate the distribution of the width of the quasienergy Wannier-Stark 
resonances and distribution of the Wigner delay time and compare them with the prediction of 
RMT. We would like to note that presently there are just a few physical models which allow a 
detailed comparison with analytical results of RMT.f] In this context, the driven Wannier-Stark 
system (|5.2D, 



H = — + cos[x + e cos(a;t)] + Fx , e = — 

^ (jj 



(7.1) 



(to be concrete, we choose V{x) = cosx) serves an excellent example for testing an abstract 
RMT. 



Among the physical m odels^ two -dimensional billiards with attached leads 1 203 -205], si mpli f ied m odels of atomic 
and molecular syst ems [206-208|, the kicked rotor with absorbing boundary condition [ 20E| -211 , and scattering 
on graphs [212,213 could be mentioned. 
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7.1. Classical dynamics 

We begin with tlie analysis of the classical dynamics of the driven Wannier-Stark system 



(7.1). Let us consider first the case F = 0. Expanding the space- and time-periodic potential in 



a Fourier series yields 

cos[x + e cos(a't)] = Jo(e) cos(x) — Ji(e) [sin(x + wt) + sin(x — a;t) ] 

— J2(e) [cos(2; + 2LJt) + cos{x — 2ujt) ] 

+J3(e) [sin(x + 3a;t) + sin(x - ^ujt)] + ... (7.2) 



Then, from the perspective of the classical nonlinear dynamics [ pl4| , the system ( ^.2[ ) is a system 
of many interacting nonlinear resonances. Depending on a particular choice of the parameters u) 



and e, its dynamics can be either quasiregular or chaotic [pl5| . This is exemplified by Fig. 7T 
where the stroboscopic surface of section^ is shown for lo = 10/6 and e = 0.1 and e = 1.5. In 
the quasiregular case with e = 0.1 only the three terms, cosx and sin(xibwt), in series ( [7.2[) are 
important. The three corresponding nonlinear resonances are clearly visible in the left panel. 
The main resonance of cosx appears as the large central island and the two other resonances 
correspond to the two smaller islands at x ~ 1 and p ~ ±1.5. For large e = 1.5 many such 
nonlinear resonances overlap, and a broad chaotic band appears. Assuming an initial condition 
in this chaotic band, the classical motion is then confined to this chaotic region, i.e. in the field 
free case F = it remains bounded in the momenta. 

Adding a dc field changes this property, since it destroys the invariant curves separating the 
chaotic component of the phase space from the outer region of the regular motion. In fact, the 
static field connects the regions of large momentum, because a particle initially localized in the 
regular region of large positive momentum p ^ p* ~ 5 can then move into chaotic region (small 
momentum \p\ < p*) from where it can finally reach the region of large negative momentum. 
Thus the scattering process p —p consists of tree stages: almost uniformly deaccelerated 
motion for p > p*, temporal chaotic motion \p\ < p*, and accelerated motion for p < —p* (see 



Fig. 7.2). The time spent by the particle in the chaotic region is the delay or dwell time r, which 



we define as the time gain or loss relative the case F(x, t) = 

r = lim [t{po ^ -po) - 2po/F] . (7.3) 

P0-+OO 



Figure 7.3 shows the delay time (measured in periods T^^) as function of the initial coordinate 
xo (the momentum pQ is kept fixed). The function is very irregular. Regions where it is approx- 
imately constant are intermitted by regions of irregular peak structures. If we zoom into such 
a structure, this behavior repeats on a finer scale, and altogether the function r(xo) shows a 
fractal behavior which is one of the main characteristics of classical chaotic scattering. 



The randomness of r suggests its statistical analysis. Figure 7A shows the distribution PdiT) 
of the classical delay time for u = 10/6, e = 1.5 and F = 0.065. It is seen that the distribution 
has an exponential tail 

Pd(^) ~ exp(-i/r) , (7.4) 

^The stroboscopic surface of section is generated by plotting the momentum p{t) and coordinate x{t), taken by 
modulus 2tt, for t = nT^ (n = 0, 1, . . .). 
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Figure 7.1. Classical stroboscopic surface of section of the system ( |7. j| j for F = 0. The driving 
frequency is uj = 10/6, and the amplitude e = 0.1 (a) and e = 1.5 (h). In the first case the 
system is almost regular, in the second case a broad chaotic band appears. 



a- Q 




Figure 7.2. Example of classical trajectories for the system ( |7. i| j with parameters oj = 10/6, 
F = 0.13 and e = (a), e = 1.5 (b). 
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Figure 7.3. Fractal structure of the classical delay time t as a function of the initial coordinate 
X. The system parameters are F = 0.3, uj = 10/6 and e = 1.5. 
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which is another characteristic feature of the chaotic scattering^. The value of the decay in- 
crement V primarily depends on F, and for F = 0.13 and F = 0.065 (used later on in the 
quantum simulation) it is ~ 0.13F or ~ 0.20F, respectively. Note that the distribution of 
the delay times also defines the decay of the classical survival probability Pci{t). Assuming an 
ensemble of classical particles with initial conditions in the chaotic region, the latter quantity 
is defined as the relative number of particles remaining in the chaotic band. Obviously, the 
classical survival probability (asymptotically) decreases exponentially with the same increment 
ly, i.e. Pci{t) PS exp(-i^t). 



7.2. Irregular quasienergy spectrum 

We proceed with the quantum mechanical analysis of the system. Let us recall that we 
consider the commensurate case of a rational ratio between the Bloch period Tb and the period 



of the exciting force, i.e. Tb/T^^ = hw = p/q with integers p and q (see Sec. |5.lD . We begin 
with the analysis of the complex quasienergy spectrum for the simplest case p = q = 1, where 
the quasienergy spectrum coincides with the spectrum of the Floquet-Bloch operator (^.^). 



^ In principle, the far asymptotic of the distribution Pci{t) may deviate from the exponential law, which is known 
to be due to the effect of the stability islands or their remnants. In our case, however, we did not observe such a 
deviation. 
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Figure 7.6. Eigenvalues Aa(«;) = exp(—i£a{K) / F) in a polar plot. The left panel shows the 
location of the eigenvalues inside the unit circle, the right panel additionally shows the dependence 
on the quasimomentum. 
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Figure 7^ shows the real and imaginary part of the spectrum <Sq,(k) for uj = 10/6, e = 1.5 
and h = 0.5. (The value of the static force is fixed by the resonant condition uj = ojb = 
I'kF jh^ which corresponds to -F ~ 0.133.) For each value of the quasimomentum the 15 most 
stable resonances are plotted. In addition to Fig. [7.5| , Fig. 7.6 shows the same spectrum as a 
polar plot for the eigenvalues Aq,(k) = exp(— iiSQ,(K)/F), where the axis of the cylinder is the 
quasimomentum axis. Now 30 resonances are depicted. It is seen in the figures that, apart from 
the symmetry k — > —k [which reflects the symmetry t — > — t, p — > — p of the Hamiltonian ( [7.1D ], 
the spectrum looks very irregular. The formal reason for this irregularity is the interaction of the 



quasienergy bands discussed in Sec. 5.4. However, in the presently considered case of small ^, 



this interaction appears to be so strong, that it makes an analytic description of the dispersion 
relation impossible. 

An important result following from the numerical data is a clear separation of the resonances 
according to their stability. Namely, for every k there is a finite number of relatively stable 



resonances which occupy the region near the unit circle in Fig. 7.6. The rest of the resonances 



are very unstable and they occupy the region in the center of the unit circle. Using the phase- 
space representation of the resonance wave function (for example the Husimi representation 
[ |161| ) it can be shown that the former resonances are supported by the chaotic region of Fig. 7T 
and, thus, are associated with the chaotic component of the classical phase space. The latter 
resonances are associated with the outer regular region of the classical phase space and can 
be considered as a kind of "above-barrier" resonances. According to the Weyl rule, the total 
number of the relatively stable (chaotic) resonances can be estimated as 



N 



1 

27r^ 



pdx , 



(7.5) 



where the integral § pdx stands for the volume of the chaotic component.^ Let us also note that 
these resonances have the width of nearly the same order of magnitude. This fact and the avoided 
crossings in the real part reflect the chaotic structure of the interaction region in classical phase 
space, which quantum mechanically results in a strong interaction of the participating states. 

Additional information about the structure of the quasienergy spectrum can be obtained by 
considering the Wigner delay time ( ^.25 ). As an example. Fig. 7.7 shows the Wigner delay time 
r(£') for u)/lob = 1/2. As already mentioned in Sec. 3.4, since 



r{E) ~ 



{E 



+ r2/4 



(7.6) 



the Wigner delay times reveals only the narrow resonances. The majority of these resonances can 



be identified with the chaotic resonances, which form an irregular pattern in Fig. 7.7. However 



besides this irregular pattern, a regular one in the form of a rhombus is clearly seen. Below we 
show that this regular structure is due to the stability islands of the classical phase space. 



In fact, let us consider an arbitrary term in Eq. (7^). This term corresponds to classical 
nonlinear resonance at p ~ muj. Assuming that the interaction between the nonlinear resonances 



"^This formula also estimates the number of under-barrier resonances for e = 0. Then j> pdx is the phase volume 
confined by the separatrix. 
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Figure 7.7. The Wigner delay time (left panel) and the real part of the quasienergy spectrum 
(right panel) for the system ( [7. i|j with parameters e = 1.5, uj = 10/6 and h = 0.25 in the case 
p/q = 1/2. 



does not completely destroy this particular resonance, the dynamics of the system in the vicinity 
of its stable periodic point is locally governed by the effective Hamiltonian 



off 



p 



+ Jm{^) cos(x lb mu!t) + Fx 



(7.7) 



(the sign of the Bessel function and the sin or cosine dependence does not matter). By sub- 
stituting x' = X i: mujt, the Hamiltonian ( [7. 71 ) is transformed to the time- independent form 
H'^g = (p^mio)'^ /2 + Jm{s) cosx' + Fx'. The latter Hamiltonian can support localized Wannier- 
Stark states^ ^ a,i{x') or, alternatively, extended Wannier-Bloch states 4?^ ,j(x'). Denoting by S'^ 
the degenerate band of these Wannier-Bloch states, the dispersion relation for the quasienergy 
spectrum of the effective Hamiltonian (|7.7|) reads (up to an additive term) 



Sain) 



[mujy 



J mod:^ 



[mio) 



+ 8'\ ±n rnlnF 



(7.8) 



mod:?ja; 



It follows from the last equation that the nonlinear resonance index m can be extracted from 
the slope of the dispersion lines. In particular, one can clearly identify the stability islands with 
m = ±2 and the remnant of the stability islands with m = ±1 in Fig. 7.1. 

To summarize, the quasienergy spectrum of the Wannier-Stark system consists of two com- 
ponents, associated with the regular and chaotic components of the classical phase space. The 
"chaotic" component of the spectrum shows a rather complicated structure. This suggests a 



'Note that these states move m lattice periods to the left or right per period of the driving frequency. 
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statistical analysis of the spectrum, which will be done in Sec. T_A. Before doing this, however, 
we shall briefly discuss some results of random matrix theory. 

7.3. Random matrix model 

As was mentioned in the introductory part of this chapter, the main conjecture of random 
matrix theory of quantum chaos is that the spectral statistics of a classically chaotic system 
coincides with those of an appropriate ensemble of the random matrix. Let us first discuss 
which ensemble is "appropriate" to model the spectral statistics of the system of our interest. 

According to the results of Sec. f>.% the quasienergy resonances of the Wannier-Stark system 
are given by the eigenvalues of the nonunitary matrix ( ^.17 ), which enters in the definition of 



the scattering matrix ( |5.23| ). In the random matrix approach it is reasonable to keep the same 
structure of the matrix. In other words, we model the case of rational u)/ojb = p/q by the 
random scattering matrix 




= eM[i?-e-'%]-V^ eM = {ONM 1m,m ) , = ( ^ ' |, (7-9) 

where nonunitary matrix B is given by 
/ Om,n Om,m \ 

B=[ ' , (7.10) 

V Wn,n On,m I 



and Wm,n is a random unitary matrix of size N x N. In Eq. ( \l.9\ ) and Eq. ( |7.10 ), the parameter 



M is identical with the denominator q in the condition of comensurability.0 Moreover, we choose 
Wn,n to be a member of Circular Unitary Ensemble (CUE). The reason for this is that matrix 
Wn,n should model the unitary matrix U^'^\ 



U^^^ = gip ( -- / dt 



T 







{p + hK- Ft) 



2 



+ V{x,t) 



(7.11) 



which, excluding the cases k = and n = ±1/2, has no time-reversal symmetry. 

Now we discuss the statistics of the resonance widths. The histograms in the left panel of 
Fig. |7.8| show the distribution of the scaled resonance widths for the random matrix model ( |7.9D , 
([7l0|)for M = 1,2,3. These histo grams are obtained in the following way. First, we generate 
a random 40 x 40 GOE (Gaussian Orthogonal Ensemble) matrix, i.e., a symmetric matrix 
with Gaussian-distributed random elements. Then, multiplying the eigenvectors of this matrix 
(arranged column-wise in a square matrix) by a random-phase factor, we obtain a member of 



CUE [216|. This CUE matrix is enlarged to a nonunitary matrix B and diagonalized. After 
diagonalization, we have (A^ — M) non-zero eigenvalues A = exp[— i(£^ — ir/2)]. To ensure the 
convergence in the limit — > oo, the resonance widths F are scaled based on the mean level 
spacing A = 27r/(A - M) as F^ = AF/2 « ttF/A. Finally, the distribution of the scaled widths 
is calculated for an ensemble of 1000 random matrices. 



In this section, we use the standard notation of RMT, i.e. TV for the matrix size and M for the number of 
scattering channels. 
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Figure 7.8. Distribution of the resonance widths (left panel) and distribution of the sum of partial 
delay time (right panel) for M = 1, 2, 3 decay channels. Numerical data (histograms) are com- 
pared with theoretical curves (solid lines). With increasing M , the maxima of the distributions 
shift to the right. 



In Fig. [7.8| the distribution of the resonance widths is compared with the analytical expression 



n(r.) 



-1 



dr 



M 



1 — ex.p{—2Ts 
2E 



(7.12) 



(M-1)! 

valid in the hmit N ^ oo Note that the distribution (|7l^ ) was originally obtained for a 

different random matrix model, which was aimed to model the chaotic scattering of the ballistic 
electrons in the mesoscopic cavities [200|, and corresponds there to the so-called case of perfect 
coupling [ [L99[ , which is realized in the case considered here. The asymptotic behavior of the 
distribution ( |7l2D is given by U{Ts) « M/2r2 for F^ > 1, and U{Ts) ~ Ff for F^ < 1. A 
perfect coincidence between the depicted numerical data and analytical results is noticed in all 
three considered cases. 

We proceed with the distribution of the Wigner delay time. The advantage of the Wigner 
delay time is that it can be directly compared to classical delay time ( |7.3| ). Within the random 
matrix approach discussed above, the Wigner delay time can be calculated by taking the trace 
of the Smith matrix (5.25), where the random matrix analogue of the scattering matrix ( ^.23| ) is 
given in Eq. ( [7. 91) . Alternatively, we can calculate the Wigner delay time by using an M-channel 
analogue of Eq. ( p.3lD 

1 



t{E) 



M 



Tr 



[5t 



e'^ll 



-lE 



1] 



(7.13) 



Note that the scattering matrix (^) yields only positive delay times whereas the Wannier-Stark 
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system (where the delay time is compared to the "free" motion) also allows for negative values. 
However, we can easily take this fact into account by shifting the delay time ( [7.13| ) by N units. 
The distributions Pqu{T) of the Wigner delay times require an additional remark. The random 



matrix theory predicts only the distribution of the partial delay times [see Eq.(7.14) below], 
whereas we are interested in the Wigner delay time, which is the sum of partial delay times 
divided by the number of channels. Because the partial delay times are correlated, the exact 
distribution of the Wigner delay time is a rather complicated problem in random matrix theory 
[218|. However, "in the first order approximation", the correlation of the partial delay times 
may be neglected. According to [219|, the correlation between partial delay times decrease as 
1/(M + 1) with increasing number of scattering channels. 
Then the distribution Pqu(T) of the Wigner delay time is the M-fold convolution of the 



distribution P{t) for the partial delay times. According to the results of Ref. |15S,199|, the 
latter is given by 



1 

M! 



(7.14) 



where Tg = t/N « rA/27r is the scaled delay time. 



In the right panel of Fig. 7.5, the distributions of the sum of the partial delay times are 



compared with the M-fold convolution of the distribution ( 7.14| )f] An ensemble contains 1000 
random matrices of the size 40 x 40, and for each matrix the delay time is calculated at 100 
equally spaced values of E. In the one channel case both results agree perfectly, whereas in the 
other cases the curves are slightly shifted. However, even here the agreement is pretty good. 
Thus, the assumption of independent partial delay times really yields a good approximation to 
the data. 



7.4. Resonance statistics 

In the previous section we introduced a random matrix model of the driven Wannier-Stark 
system which yields analytical results for the distribution of the resonance width and Wigner 
delay time. In this section we compare the actual distributions, obtained numerically, to these 
theoretical predictions. In our calculation, we construct the statistical ensemble by scanning the 
quasimomentum k with a step Ak over the first Brillouin zone —l/2p < k < l/2p. To get a 
good statistics, Ak should be as small as possible. On the other hand, because the widths and 
the delay times depend smoothly on the quasimomentum, there is a characteristic value of Ak 
such that a further decrease does not improve the statistics. In the following calculations we 
choose Ak = l/200p, i.e. we average over 200 spectra. 

The other problem arising in the statistical analysis of the numerical data is the appropriate 
rescaling of the resonance width and the delay time. In fact, the notion of matrix size A'^ is 
not directly specified in our approach. However, we can use the semiclassical estimate (|7.5| ) to 
specify the parameter N. For the value of the scaled Plank constant h = 0.25 considered below 
this gives ~ 32. In what follows, however, we use a slightly smaller value A^ = 28 which 
accounts for the embedded islands of stability. 



'^To obtain the distribution for the Wigner delay time, the displayed histograms should be scaled as Pqu(r) 
MP{t/M). 
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Figure 7.9. Complex quasienergies, distribution of the widths, delay time and distribution of the 



delay times of the system (7.1) for the case p/q = 1 with parameters e = 1.5, to = 10/6 and 
h = 0.25. In this case the constant force is F 0.066. 



The two upper panels of Fig. |7.9| show the complex quasienergies and the distribution of the 
resonance widths for the most simple single channel case p = q = 1. A good agreement between 
the random matrix results and the calculated distribution is noticed. The distribution has its 
maximum at F^ = 0, i.e. the resonances tend to be long-lived. The main deviation is a peak at 
Tg ~ 8, which is due to states associated with stability islands in the classical surface of section. 



One also finds these resonances in the delay time shown in the left lower panel of Fig. 7.9. As 



discussed in Sec. 7.2, resonances corresponding to classical stability islands form straight lines 
in the quasienergy spectrum. Indeed, we can see the lines with the slope ib47r and remnants of 
two lines with slope ib27r. Because such resonances have approximately the same widths, their 
signatures are easily identified in the distribution of the widths. 

The right lower panel shows the distribution of the delay time (to facilitate the comparison, 
the histogram for the scaled delay time is shifted to the right by one unit.) Here the agreement is 
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Figure 7.10. Distribution of the resonance widths for different field strengths F with e = 1.5, 
uj = 10/6 and h = 0.25. The histogram show the numerical data, the solid lines are the random 
matrix predictions { 7.1i ) for the proper number of decay channels. 
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pretty good, either. The location of the maxima at = 0.33 and the shape of both distributions 
coincide almost perfectly. 

We proceed with the case p/q ^ 1, where we restrict ourselves to an analysis of the resonance 
widths. The most striking prediction of the random matrix model of Sec. |7.3| is that the statistics 
of the resonance widths is solely defined by the integer q. On the other hand, the random matrix 
model is supposed to describe the properties of the real system with four parameters. Thus, 
provided q is the same, the distribution of the resonance widths should be independent on the 
particular choice of the other system parameters. (Of course, the condition for chaotic dynamics 
should be fulfilled.) To check this prediction we proceed as follows. 

The number q of decay channels is defined by the rationality condition pT^^ = qTs, i.e. in 
terms of the system parameters hy F = qhjj/p2Ti. As in the preceding cases, we choose e = 1.5 
and ljo = 10/6 to ensure that the system is classically chaotic, and h = 0.25 in order to be in 
the semiclassical regime. Then we calculate the distribution of the resonance widths for several 
combinations of the integers p and q, which correspond to increasing values of F. Naively one 
would expect that with increasing F the resonances tend to destabilize. Instead the distributions 
follow closely the RMT distributions for the g-channel case as can be seen in Fig. 7.1C. The first 
and the last picture for the smallest and the largest field strength correspond to the one-channel 
case q = 1. Note that the field strength differs by a factor seven, but the distributions are 
essentially the same. For the intermediate field strengths the distributions vary according to the 
number of decay channels. We should stress that the only adjusted parameter, the number of 
states N = 28 defining the scaled width r<j, is constant in all figures. 



7.5. Fractional stabilization 

In this section we discuss an interesting application of the results of preceding section, which 
can be referred to as fractional stabilization of the Wannier-Stark system. 

Let us discuss again the spectroscopic experiment |123|] , where the survival probability for the 
cold atoms in the accelerated optical lattice was measured as a function of the driving frequency 
(see Sec. [1.3|).We assume now the following modifications of the experimental setup. The value 
of the scaled Planck constant (which is inversely proportional to the laser intensity) is small 
enough to insure the semiclassical dynamics of the system. The value of the driving amplitude 
is large enough to guarantee the classical chaotic dynamics of the atoms. (Note that both these 
condition were satisfied in a different experiment p21[| .) The atomic survival probability is 
measured as a function of the acceleration but not as the function of the driving frequency, i.e. 
we vary ujb instead of varying u). (This condition is actually optional.) 

Figure [7.11 shows the results of the numerical calculation of the quantum survival probabil- 



ity Pqu{t) based on direct numerical simulation of the wave packet dynamics. ^ The survival 
probability shows an interesting behavior. For small times, the curve fiuctuates around an 
approximately constant value. When the time is increased, this average value decreases expo- 
nentially. In addition, however, peaks develop at integer values 7 = uj/ujb and, incrementally, 



^Explicitly, we calculate wave function 'ip{p,t) in the momentum representation with the localized Wannier state 
as an initial condition. Then the probability for a quantum particle to stay within the chaotic region is given by 
Pqu{t) — i|p|<pt IV'(PiOPdp, where p* is the classical boundary between the chaotic and regular components of 
the classical phase space. 
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Figure 7.11. Survival probability as function of the control parameter 7 = lj/lob- The system 
parameters are uj = 10/6, e = 1.5 and h = 0.25. 



at rational "f = p/q with small denominator. Thus, the decay is slowed down for rational 7. In 
what follows we explain this stabilization effect by using RMT approach. 

Indeed, the system parameters were chosen to ensure the regime of chaotic scattering. Then 
the distribution of the resonance widths is given by equation ( [7.12| ). Let us assume that the 
initial state uniformly populates all resonances. If we then neglect the overlap of the resonances 
(this is the so-called diagonal approximation) the survival probability is given by the integral 

/•oo 

p^.{t)= / drn(r)e-r*/\ (7.15) 

Jo 



where T = 2Ts/N and N is the number of states in the interaction region. The long-time 



82 



aoB 



008 




Figure 7.12. Survival probability at t = 200Tb. The left part shows the numerical data, the right 
part the theoretical curve based on equation (7.16). To stress the discontinuous character of the 
latter function, we slightly changed its graphical representation. 



asymptotics of this integral is defined by the behavior of n(r) at small F, where it increases as 
the power law n(F) ~ T'^~^. Consequently, the survival probability asymptotically follows the 
inverse power law Pqu{t) ~ i Thus, the asymptotics depend on the number of decay channels 
and therefore on the denominator of the control parameter 7 = uj/ujb = p/q- 

With the help of supersymmetric techniques, Pquit) can be calculated beyond the diagonal 
approximation. This gives more elaborated result ]22C| ] 

where Tw is the so-called Weisskopf width (which is a free parameter in the abstract ran- 
dom matrix theory). For rational 7 and large times the decay of the survival probability is 
algebraic, -Pqu(i) ~ t^^, as found in the diagonal approximation. The case of irrational 7 
can be approximated by the limit q ^ 00. Then the system shows the exponential decay, 
-fqu(i) = exp(— Fiyt/^) and its natural to identify the parameter Tw/fi with the classical decay 
coefficient v. 

The right panel in Fig. 7.12 shows the values of the function ( [7.16| ) for t = 200 Tg and some 
rational values of 7 = uj/lob- Here we use a slightly different graphic presentation of Pquit) to 



stress that the function (7.16) is a discontinuous function of 7 for any t. In contrast, the atomic 



survival probability shown in the left panel is a continuous function of 7 where its discontinuous 
structure develops gradually as t — > 00. In fact, the survival probabilities calculated for two 
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close rational numbers 71 and 72 follow each other during a finite "correspondence" time. (For 
instance, for 71 = 1 and 72 = 999/1000 the correspondence time is found to be about 50 Tg.) 
Thus it takes some time to distinguish two close rationals, although they may have very different 
denominators and, therefore, very different asymptotics. With this remark reserved, a nice 
structural (and even semiquantitative) correspondence is noticed. 

The described numerical experiment suggest a simple laboratory experiment with cold atoms 
in optical lattice, where one can test the statistics of the resonance width indirectly, by measuring 
the survival probability for atoms. 



Chapter 8 

Conclusions and outlook 



In this section we review the main results of the work and outline some problems which are still 
waiting for their solutions. In the overview we shall mainly follow the table of contents. 

The approach introduced in Sec. 2 gives us a powerful tool for analyzing an arbitrary one- 
dimensional Wannicr-Stark system, i.e. a system with potential energy given by the sum of 
periodic and linear terms. The success of the method is ensured by two key points. First, 
we inverted the traditional solid state approach, where the linear term has been treated as 
"a perturbation" to the periodic term, and formulate the problem as scattering of a quantum 
particle by a periodic potential. Second, instead of dealing with the Hamiltonian, we work 
with the evolution operator. Although both these points were discussed earlier, it is only a 
combination of them, which provides solution of the Wannier-Stark problem. Let us also note 
that the use of the evolution operator provides a way to an analysis of the Wannier-Stark 
system affected additionally by a time-periodic perturbation. The corresponding generalization 
of the method, which leads to the notion of the metastable quasienergy Wannier-Stark states, 
is discussed in Sec. 5 of this review. 

We apply the developed theory to analyze the Wannier-Stark ladder of resonances in two 
particular systems - undoped semiconductor superlatticcs in a static electric field and the system 
of cold atoms in optical lattices in an accelerated frame. Both of these systems mimic the crystal 
electron in a static electric field (which was the original formulation of the problem) and have 
their own advantages and disadvantages. In particular, the semiconductor superlattices allow 
(at least, in principle) to create an arbitrary periodic potential. One may think, for example, 
about a periodic sequence of double wells, where the interaction of the Wannier-Stark ladders 
(which is essentially the resonance tunneling effect) should have an especially interesting form. 
In Sec. 3 we restricted ourselves by considering the cosine and square-box shaped potentials. The 
structure of the Wannier-Stark states and the interaction of the ladders in periodic potentials 
of a different form (like the already mentioned double-well array or asymmetric ratchet-like 
potential) is an open problem. 

A disadvantage of the semiconductor superlattice is that this is a more "dirty" (in comparison 
with the optical superlattice) system, where the effects in question interfere with other effects 
like electron-hole Coulomb interaction, scattering by impurities, etc. Nevertheless, if we want 
to move further, we should learn how to deal with these complications. In the first turn, the 
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effect of Coulomb interaction should be taken into account. We believe that now this problem 
can be solved rigorously by extending the one-particle scattering theory of Sec. 2 to the case of 
two particles. 

We turn to the spectroscopic results of Sec. 4. In this section we derive an analytic expression 
for the decay spectrum of the system of cold atoms in an optical lattice and the absorption spec- 
trum of semiconductor super lattices. This expression involves complex valued squared transition 
matrix elements (non-real squared matrix elements appear because of the resonance nature of 
the Wannier-Stark states), which lead to a non-Lorentzian shape of the absorption lines. Al- 
though the relation of these results to the famous Fano theory is obvious, the details of this 
relation remain unexplored. 

The brief Sec. 6 was inspired by the experiment of Anderson and Kasevich, where a pulsed 
output from the periodic array of cold atoms was observed. We give a proper theoretical de- 
scription of this phenomenon which, in fact, is the Bloch oscillations in the case of a strong 
static field. In this sense. Sec. 6 is the only section of the review discussing Bloch oscillations. 
One might be interested in other regimes of Bloch oscillations. Evolution of the theory in this 
direction is reflected by a recent paper |222| |. 

As already mentioned in the introduction. Sec. 7 deals with the very different problem of 
chaotic scattering, which is primary of interest to the members of quantum chaos community. 
Nevertheless, from the formal point of view, the results of Sec. 7 are just the results beyond the 
perturbative approach of Sec. 5. Thus, when the experimentalists overcome the perturbation 
limit (the present state of the art), Sec. 7 may change its status from of "pure theoretical 
interest" to that of "practical importance". 

To conclude, we would like to highlight one more problem. This work is devoted entirely 
to one-dimensional Wannier-Stark systems. However, practically nothing is known about the 
Wannier-Stark states in 3D- or 2D-lattices (a first step in this direction was taken only recently 
[p23| ). An extension of the present theory to higher dimension is of much theoretical and 
practical interest and one may expect on this way a variety of new phenomena which are absent 
in the one-dimensional case. 

Furthermore, the results presented in this review will also be relevant in connection with 
recent new developments in quantum transport in driven periodic lattices with broken symmetry, 
i.e. quantum hamiltonian ratchets |224-228|. Such ratchets are usually studied in the case of 
vanishing mean potential gradient. An interesting situation arises, e.g., for ratchets inclined 
in the direction opposed to the current that would occur in the unbiased case. In addition, it 
should be noted that in the previous studies of the classical - quantum correspondence for driven 
Wannier-Stark systems as discussed in Sec. |^ the parameters have been chosen to guarantee 
(almost) fully chaotic dynamics in the scattering region, i.e. classical stability islands are of 
minor importance. Larger islands can be observed, however, and will certainly effect the decay 
properties discussed in Se c. ^ , as for instance by chaos-assisted tunneling, a topic of much 
interest in theoretical [229-231] and very recently also experimental studies 1 232-234], were this 
phenomenon was rediscovered. 



86 



Acknowledgments 

The authors gratefully acknowledge the support of the Deutsche Forschungsgemeinschaft via 
the Schwcrpunktprogramm SSP 470 "Zeitabhangigc Phanomene und Methoden in Quanten- 
systemen der Physik und Chemie" as well as the Graduiertenkolleg "Laser- und Teilchenspek- 
troskopie". Stimulating discussions with many colleagues are gratefully acknowledged, in par- 
ticular, among others, with J. E. Avron, Y. V. Fyodorov, M. Holthaus, K. Leo, N. Moiseyev, 
Q. Niu, M. G. Raizen, D. V. Savin, V. V. Sokolov, and K. Zyczkowski. It is also a pleasure to 
thank the graduate students Michael Hankel, Christian Hebell, Frank Keck, Stefan Mossmann, 
Prank Zimmer who contributed to our studies of Wannier- Stark systems. Moreover, we thank 
P. Hanggi for many useful suggestions which considerably improved the present review. 



87 



REFERENCES 

1. F. Bloch, Uber die Quantenmechanik der Electronen in Kristallgittem, Z. Phys 52 (1928) 

555 

2. C. Zener, A Theory of Electrical Breakdown of Solid Dielectrics, Roy. Soc. Lond. A 145 
(1934) 523 

3. G. H. Wannier, Wave Functions and Effective Hamiltonian for Bloch Electrons in an Electric 
Field, Phys. Rev. 117 (1960) 432 

4. J. Zak, Stark Ladder in Solids?, Phys. Rev. Lett. 20 (1968) 1477 

5. G. H. Wannier, Stark Ladder in Solids? A Reply, Phys. Rev. 181 (1969) 1364 

6. J. Zak, Stark Ladder in Solids? A Reply to a Reply, Phys. Rev. 181 (1969) 1366 

7. A. Rabinovitch and J. Zak, Electrons in Crystals in Finite-Range Electric Fields, Phys. Rev. 
B 4 (1971) 2358 

8. W. Shockley, Stark Ladders for Finite, One-Dimensional Models of Crystals, Phys. Rev. 
Lett. 28 (1972) 349 

9. A. Rabinowitch and J. Zak, Does a Bloch Electron in a Constant Electric Field Oscillate?, 
Phys. Lett. A 40 (1972) 189 

10. J. N. Churchill and F. E. Holmstrom, Comments on the Existence of Bloch Oscillations, 
Phys. Lett. A 85 (1981) 453 

11. J. N. Churchill and F. E. Holmstrom, Energy States and Bloch States for an Accelerated 
Electron in a Periodic Lattice, Phys. Scripta 27 (1983) 91 

12. J. B. Krieger and G. J. lafrate, Time Evolution of Bloch Electrons in a Homogeneous Electric 
Field, Phys. Rev. B 33 (1986) 5494 

13. D. Emin and C. F. Hart, Existence of Wannier-Stark Localization, Phys. Rev. B 36 (1987) 
7353 

14. C. F. Hart and D. Emin, Time evolution of a Bloch electron in a constant electric field, 
Phys. Rev. B 37 (1988) 6100 

15. L. Kleinman, Comment on "Existence of Wannier-Stark localization", Phys. Rev. B 41 

(1990) 3857 

16. J. Zak, Comment on the existence proofs of the Wannier-Stark ladder, Phys. Rev. B 43 

(1991) 4519 

17. D. A. Page and E. Brown, Comment on "Existence of Wannier-Stark localization", Phys. 
Rev. B 43 (1991) 2423 

18. J. Leo and A. MacKinnon, Comment on "Existence of Wannier-Stark localization", Phys. 
Rev. B 43 (1991) 5166 

19. X. G. Zhao, Bloch Electron in a Uniform Electric Field, Phys. Rev. B 46 (1992) 1305 

20. G. Nenciu, Dynamics of hand electrons in electric and magnetic fields: rigorous justification 
of the effective Hamiltonians, Rev. Mod. Phys. 63 (1991) 91 

21. A. M. Bouchard and M. Luban, Bloch Oscillations and other Dynamical Phenomena of 
Electrons in Semiconductor Superlattices, Phys. Rev. B 52 (1995) 5105 

22. F. Rossi, Bloch Oscillations and Wannier-Stark Localization in Semiconductor Superlattices, 
in E. Scholl, editor, Theory of Transport Properties of Semiconductor Nanostructures, page 
283. Chapman & Hall, London, 1998 

23. J. E. Avron, J. Zak, A. Grossmann, and L. Gunther, Instability of the continuous spectrum: 



88 



The N-band Stark ladder, J. Math. Phys. 18 (1977) 918 

24. F. Bentosela, R. Carmona, P. Duclos, B. Simon, B. Souillard, and R. Weder, Schrddinger 
Operators with an Electric Field and Random or Deterministic Potentials, Com. Math. Phys. 
88 (1983) 387 

25. I. W. Herbst and J. S. Rowland, The Stark Ladder and Other One-Dimensional External 
Field Problems, Com. Math. Phys. 80 (1981) 23 

26. J. Agler and R. Proese, Existence of Stark Ladder Resonances, Com. Math. Phys. 100 (1985) 
161 

27. J.-M. Combes and P. D. Hislop, Stark Ladder Resonances for Small Electric Fields, Com. 
Math. Phys. 140 (1991) 291 

28. F. Bentosela and V. Grecchi, Stark Wannier Ladders, Com. Math. Phys. 142 (1991) 169 

29. N. Moiseyev, Quantum, theory of resonances: calculating energies, widths and cross-sections 
by complex scaling, Phys. Rep. 302 (1998) 211 

30. J. E. Avron, Model Calculations of Stark Ladder Resonances, Phys. Rev. Lett. 37 (1976) 
1568 

31. J. E. Avron, The Lifetime of Wannier Ladder States, Ann. Phys. (N.Y.) 143 (1982) 33 

32. V. Grecchi, M. Maioli, and A Sacchetti, Wannier Ladders and Perturbation Theory, J. Phys. 
A 26 (1993) L379 

33. V. Grecchi, M. Maioh, and A Sacchetti, Stark Ladder of Resonances: Wannier Ladders and 
Perturbation Theory, Com. Math. Phys. 159 (1994) 605 

34. V. Grecchi and A. Sacchetti, Crossing and Anticrossing of Resonances: The Wannier-Stark 
Ladders, Ann. Phys. (N.Y.) 241 (1995) 258 

35. V. Grecchi and A. Sacchetti, Metastable Bloch Oscillators, Phys. Rev. Lett. 78 (1997) 4474 

36. V. Grecchi and A. Sacchetti, Lifetime of the Wannier-Stark Resonances and Parturbation 
Theory, Com. Math. Phys. 185 (1997) 359 

37. V. Grecchi and A. Sacchetti, Wannier-Bloch Oscillators, Com. Math. Phys. 197 (1998) 553 

38. V. Buslaev and A. Grigis, Lmaginary Parts of Stark- Wannier Resonances, J. Math. Phys. 
39 (1998) 2520 

39. J. Banavar and D. D. Coon, Widths and Spacing of Stark Ladder Levels, Phys. Rev. B 17 
(1978) 3744 

40. F. Bentosela, V. Grecchi, and F. Zironi, Approximate ladder of resonances in a semi-infinite 
crystal, J. Phys. C: Solid State Phys. 15 (1982) 7119 

41. F. Bentosela, V. Grecchi, and F. Zironi, Oszillations of Wannier Resonances, Phys. Rev. 
Lett. 50 (1983) 84 

42. M. Ritze, N. J. M. Horing, and R. Enderlein, Density of States and Wannier-Stark Levels 
of Superlattices in an Electric Field, Phys. Rev. B 47 (1993) 10437 

43. M. C. Chang and Q. Niu, Local density of states and level width for Wannier-Stark ladders, 
Phys. Rev. B 48 (1993) 2215 

44. C. L. Roy and P. K. Mahapatra, Bloch electrons in finite crystals in the presence of a uniform 
electric field, Phys. Rev. B 25 (1982) 1046 

45. B. Soucail, R. Ferreira, G. Bastard, and P. Voisin, Instability of Energy Band Structure at 
Low Electric Field: Numerical Analysis of Superlattice Minibands, Europhys. Lett. 15 (1991) 
857 

46. B. Mendez and F. Dominguez-Adame, Stark Ladders in Periodically Si-S-Doped GaAs, Phys. 



89 



Rev. B 49 (1994) 11471 

47. S. Glutsch and F. Bechstedt, Interaction of Wannier- Stark Ladders and Electrical Breakdown 
in Superlattices, Phys. Rev. B 60 (1999) 16584 

48. H. Fukuyama, R. A. Bari, and H. C. Fogedby, Tightly Bound Electrons in a Uniform Electric 
Field, Phys. Rev. B 8 (1973) 5579 

49. D. H. Dunlap and V. M. Kenkre, Dynamic localization of a charges particle moving under 
the influence of an electric field, Phys. Rev. B 34 (1986) 3625 

50. X.-G. Zhao, Dynamic Localization Conditions of a Charged Particle in a dc-ac Electric Field, 
Phys. Lett. A 155 (1991) 299 

51. X.-G. Zhao, Motion of Bloch Electrons in Time- Dependent Electric Fields with Off-Diagonal 
Effects, Phys. Lett. A 167 (1992) 291 

52. N. Hhong-Shon and H. N. Nazareno, J. Phys. Cond. Matter 4 (1992) L611 

53. X.-G. Zhao and Q. Niu, Localization of Band Electrons in dc-ac Electric Fields, Phys. Lett. 
A 191 (1994) 181 

54. X.-G. Zhao, R. Jahnke, and Q. Niu, Dynamic Fractional Stark Ladders in dc-ac Fields, 
Phys. Lett. A 202 (1995) 297 

55. K. Drese and M. Holthaus, Exploring a Metal- Insulator Transition with Ultracold Atoms in 
Standing Light Waves?, Phys. Rev. Lett. 78 (1997) 2932 

56. M. Grifoni and P. H. Hanngi, Driven Quantum Tunneling, Phys. Rep. 304 (1998) 229 

57. D. W. Hone and M. Holthaus, Locally Disordered Lattices in Strong Ac Electric Fields, Phys. 
Rev. B 48 (1993) 15123 

58. H. Yamada, K. Ikeda, and M. Goda, Quantum diffusion in a coherently time-varying one- 
dimensional disordered system, Phys. Lett. A 182 (1993) 77 

59. M. Holthaus, G. H. Ristow, and D. W. Hone, ac-Field-Controlled Anderson Localization in 
Disordered Semiconductor Superlattices, Phys. Rev. Lett. 75 (1995) 3914 

60. M. Holthaus, G. H. Ristow, and D. W. Hone, Random lattices in combined ac and dc electric 
fields: Anderson vs. Wannier-Stark localization, Europhys. Lett. 32 (1995) 241 

61. K. Drese and M. Holthaus, Anderson Localization in an ac-Driven Two-Band Model, J. 
Phys. Condens. Matter 8 (1996) 1193 

62. H. Yamada, K. S. Ikeda, Anomalous diffusion and scaling behavior of dynamically perturbed 
one- dimensional disordered quantum systems, Phys. Lett. A 248 (1993) 179 

63. D. Suqing and X.-G. Zhao, Effect of External Noise on Dynamic Localization of a Charged 
Particle, Phys. Rev. B 61 (2000) 5442 

64. V. I. Kovanis and V. M. Kenkre, Exact Self- Propagators for Quasiparticle Motion on a Chain 
with Alternating Site Energies or Intersite Interactions, Phys. Lett. A 130 (1988) 147 

65. X.-G. Zhao, Exact Solutions for a Charged Particle in a Uniform Electric Field with Alter- 
nating Site Energies: Perturbation Theory, J. Phys.: Condens. Matter 3 (1991) 6021 

66. X.-G. Zhao, Dynamics of a Periodic Binary Sequence in an AC Field, J. Phys.: Condens. 
Matter 9 (1997) L385 

67. S.-Q. Bao, X.-G. Zhao, X.-W. Zhang, and W.-X. Yan, Dynamics of a Charged Particle on 

a Periodic Binary Sequence in an External Field, Phys. Lett. A 240 (1998) 7771 

68. P. H. Rivera and P. A. Schulz, Tuning of Dynamic Localization in Coupled Minibands: 
Signatures of a Field-Induced Metal-Insulator Transition, Phys. Rev. B 61 (2000) R7865 

69. J. Rotvig, A.-P. Jauho, and H. Smith, Bloch Oscillations, Zener Tunneling and Wannier- 



90 



Stark Ladders in the Time Domain, Phys. Rev. Lett. 74 (1995) 1831 

70. J. Rotvig, A.-P. Jauho, and H. Smith, Theory of coherent time- dependent transport in one- 
dimensional multiband semiconductor superlattices, Phys. Rev. B 54 (1996) 17691 

71. D. W. Hone and X.-G. Zhao, Time-periodic behaviour of multiband superlattices in static 
electric fields, Phys. Rev. B 53 (1996) 4834 

72. X.-G. Zhao, W.-X. Yan, and D. W. Hone, Zener transitions between dissipative Bloch bands, 
Phys. Rev. B 57 (1998) 9849 

73. W.X. Yan, X.-G. Zhao, and S.-Q. Bao, Dynamics of two-band semiconductor superlattices 
driven by static and time- dependent fields, Physica B 252 (1998) 63 

74. X.-G. Zhao, G. A. Georgakis, and Q. Niu, Rabi oscillations between Bloch bands, Phys. Rev. 
B 54 (1996) R5235 

75. W. Kohn, Construction of Wannier Functions and Applications to Energy Bands, Phys. 
Rev. B 7 (1972) 4388 

76. W. Kohn, Analytic Properties of Bloch Waves and Wannier Functions, Phys. Rev. 115 
(1959) 809 

77. G. Nenciu, Existence of the Exponentially Localized Wannier Functions, Com. Math. Phys. 
91 (1983) 81 

78. W. V. Houston, Acceleration of Electrons in a Crystal lattice, Phys. Rev. 57 (1940) 184 

79. L. D. Landau, Phys. Z. Sov. 1 (1932) 46 

80. R. W. Koss and L. M. Lambert, Experimental Observation of Wannier Levels in Semi- 
Insulating Gallium Arsenide, Phys. Rev. B 5 (1972) 1479 

81. L. Esaki, The Evolution of Semiconductor Quantum Structures. Do-It-Yourself Quantum 
Mechanics,, in A. Stella and L. Migho, editors, Superreticoli e Interface die Semiconduttori, 
number 117 in Proceedings of the International School of Physics "Enrico Fermi", page 1, 
Amsterdam, 1993. North-Holland 

82. E. E. Mendez, F. Agullo-Rueda, and J. M. Hong, Stark Localizations in GaAs-GaAlAs 
Superlattices under an Electric field, Phys. Rev. Lett. 60 (1988) 2426 

83. P. Voisin, J. Bleuse, C. Bouche, S. Gaillard, C. Alibert, and A. Regreny, Observation of 
the Wannier-Stark quantization in a semiconductor superlattice, Phys. Rev. Lett. 61 (1988) 
1639 

84. E. E. Mendez and G. Bastard, Wannier-Stark Ladders and Bloch Oscillations in Superlat- 
tices, Physics Today 46(6) (1993) 34 

85. J. Feldmann, K. Leo, J. Shah, B. A. B. Miller, J. E. Cunningham, T. Meier, G. von Plessen, 
A. Schulze, P. Thomas, and S. Schmitt-Rink, Optical Investigation of Bloch Oscillations in 
a Semiconductor Superlattice, Phys. Rev. B 46 (1992) 7252 

86. K. Leo, P. H. Bolivar, F. Briiggemann, R. Schwedler, and K. Kohler, Observation of Bloch 
Oscillations in a Semiconductor Superlattice, Solid State Comm. 84 (1992) 943 

87. G. von Plessen and P. Thomas, Method for Observing Bloch Oscillations in the Time Do- 
main, Phys. Rev. B 45 (1992) 9185 

88. C. Waschke, H. G. Roskos, R. Schwedler, K. Leo, H. Kurz, and K. Kohler, Coherent 
Submillimeter-Wave Emission from Bloch Oscillations in a Semiconductor Superlattice, 
Phys. Rev. Lett. 70 (1993) 3319 

89. T. Dekorsy, P. Leisching, K. Kohler, and H. Kurz, Electro-Optic Detection of Bloch Oscilla- 
tions, Phys. Rev. B 50 (1994) 8106 



91 



90. G. C. Cho, T. Dekorsy, H. J. Bakkcr, H. Kurz, A. Kohl, and B. Opitz, Block- Oscillations in 
In-Ga-As-P/In-Ga-As-P Heterostructures Oberved with Time-Resolved Transmission Spec- 
troscopy, Phys. Rev. B 54 (1996) 4420 

91. K. Leo, Interband Optical Investigation of Block Oscillations in Semiconductor Superlattices, 
Semicond. Sci. TechnoL 13 (1998) 249 

92. T. Dekorsy, R. Ott, and H. Kurz, Block Oscillations at Room Temperature, Phys. Rev. B 
51 (1995) 17275 

93. V. G. Lyssenko, G. Valusis, F. Loser, T. Hasche, and K. Leo, Direct Measurement of tke 
Spatial Displacement of Block-Oscillating Electrons In Semiconductor Superlattices, Phys. 
Rev. Lett. 79 (1997) 301 

94. M. Sudzius, V. G. Lyssenko, F. Loser, K. Leo, M. M. Dignam, and K. Kohler, Optical Control 
of Block-Oscillation Amplitudes: Prom Harmonic Spatial Motion to Breatking Modes, Phys. 
Rev. B 57 (1998) R12693 

95. B. Soucail, N. Dupuis, R. Ferreira, P. Voisin, A. P. Roth, D. Morris, K. Gibb, and C. Lacehe, 
Electron Minihands and Wannier- Stark Quantization in an Ino^i^Gao^sb^s-GaAs Strained- 
Layer Superlattice, Phys. Rev. B 41 (1990) 8568 

96. D. M. Whittaker, M. S. Skolnick, G. W. Smith, and C. R. Whitchousc, Wannier-Stark 
Localization of X and T States in GaAs-AlAs Skort-Period Superlattices, Phys. Rev. B 42 
(1990) 3591 

97. H. Schneider, K. Jawashima, and K. Fujiwara, Stark Localization of a Pair of Coupled 
Minihands in GaAs/AlAs Double-Period Superlattice, Phys. Rev. B 44 (1991) 5943 

98. R. H. Yu, Wannier-Stark Localization in Modulation-Doped Multiple-Quantum-Well Struc- 
tures, Phys. Rev. B 49 (1994) 4673 

99. D. W. Peggs, M. S. Skolnik, D. M. Whittaker, R. A. Hogg, A. R. K. Willcox, D. J. Mowbray, 
R. Grey, G. J. Rees, L. Hart, M. Hopkinson, G. Hill, and M. A. Pate, Observation of 
Wannier-Stark Ladder Transitions in InxGai-xAs-GaAs Piezoelectric Superlattices, Phys. 
Rev. B 52 (1995) R14340 

lOOK. Gibb, M. M. Digman, J. E. Sipe, and A. P. Roth, Observation of Wannier-Stark Local- 
ization by Electroreflectance Spectroscopy, Phys. Rev. B 48 (1993) 8156 

lOl.C. Hamaguchi, M. Yamaguchi, M. Morifuji, H. Kubo, K. Taniguchi, C. Gmachl, and 
E. Gornik, Wannier-Stark Effect in Superlattices, Semicond. Sci. Technol. 9 (1994) 1994 

102H. Schneider, H. T. Grahn, K. v. Klitzing, and K. Ploog, Resonance- Induced Delocalization 
of Electrons in CaAs-AlAs Superlattices, Phys. Rev. Lett. 65 (1990) 2720 

103M. Nakayama, I. Tanaka, H. Nishimura, K. Kawashima, and K. Fujiwara, Electroreflectance 
Detection of Resonant Coupling between Wannier-Stark Localization States in Semiconduc- 
tor Superlattices, Phys. Rev. B 44 (1991) 5935 

104.1. Tanaka, M. Nakayama, H. Nishimure, K. Kawashime, and K. Fujiwara, Electroreflectance 
Intensity for Resonant Coupling between Wannier-Stark Localisation States in GaAs/AlAs 
Superlattice, Phys. Rev. B 46 (1992) 7656 

105. G. Bastard, R. Ferreira, S. Chelles, and P. Voisin, Interaction between Wannier- States in 
Semiconductor Superlattices, Phys. Rev. B 50 (1994) 4445 

106. H. Kiimmel, R. Till, and A. Philip, Photocurrent Spectroscopy of Low- Electric- Field Anti- 
crossings in Semiconductor Superlattices, Phys. Rev. B 60 (1999) 4470 

107. G. von Plessen, T. Meier, J. Feldmann, E. O. Gobel, P. Thomas, K. W. Goossen, J. M. 



92 



Kuo, and R. F. Kopf, Influence of Scattering on the Formation of Wannier- Stark Ladders 
and Block- Oscillations in Semiconductor Superlattices, Phys. Rev. B 49 (1994) 14058 

108. J.-B. Xia, Scattering Rates of Wannier States in Superlattices in an Electric Fields Phys. 
Rev. B 50 (1994) 15067 

109. F. Loser, Y. A. Kosevich, K. Kohler, and K. Leo, Dynamics of Block Oscillation under tke 
Influence of Scattering and Cokerent Plasmon Coupling, Phys. Rev. B 61 (2000) R13373 

110. E. Ribeiro, F. Cerdeira, and A. P. Roth, Step-by-Step Evolution from Franz-Keldysh Oscil- 
lations to Wannier- Stark confinement in an Ino_i2Gao_s8 /Go-As Superlattice, Phys. Rev. B 
46 (1992) 12542 

111. K. H. Schmidt, N. Linder, G. H. Dohler, H. T. Grahn, K. Ploog, and H. Schneider, Coex- 
istence of Wannier- Stark Transitions and Minihand Franz-Keldysh Oscillations in Strongly 
Coupled GaAs-AlAs Superlattices, Phys. Rev. Lett. 72 (1994) 2769 

112. N. Linder, K. H. Schmidt, W. Geissclbrccht, G. H. D61ilcr, H. T. Grahn, K. Ploog, and 
H. Schneider, Coexistence of tke Franz-Keldysk and Wannier-Stark Effect in Semiconductor 
Superlattices, Phys. Rev. B 52 (1995) 17352 

113. M. M. Dignam and J. E. Sipe, Exciton Stark Ladders in Semiconductor Superlattices, Phys. 
Rev. B 43 (1991) 4097 

114. A. M. Fox, C. A. B. Miher, J. E. Cunningham, W. Y. Jan, C. Y. P. Chao, and S. L. Chuang, 
Suppression of the Observation of Stark Ladders in Optical Measurements on Superlattices 
by Excitonix Effects, Phys. Rev. B 46 (1992) 15365 

115. P. Leisching, P. H. Bohvar, W. Beck, Y. Dhaibi, F. Briiggemann, R. Schwedler, H. Kurz, 
K. Leo, and K. Kohler, Block Oscillations of Excitonic Wave Packets in Semiconductor 
Superlattices, Phys. Rev. B 50 (1994) 14389 

116JVI. Dignam, J. E. Sipe, and J. Shah, Cokerent Excitations in tke Stark Ladder: Excitonic 

Block Oscillations, Phys. Rev. B 49 (1994) 10502 
117JSr. Linder, Excitons in Superlattices: Absorption Asymmetry, Dimensionality Transition and 

Exciton Localization, Phys. Rev. B 55 (1997) 13664 
118 A. di Carlo, P. Vogl, and W. Potz, Theory of Zener Tunneling and Wannier-Stark states in 

Semiconductors, Phys. Rev. B 50 (1994) 8358 
119.A. Sibille, J. F. Palmier, and F. Laruelle, Zener Interminiband Resonant Breakdown in 

Superlattices, Phys. Rev. Lett. 80 (1998) 4506 
120M. Helm, W. Hilber, G. Strasser, R. de Meester, F. M. Peeters, and A. Wacker, Continuum 

Wannier-Stark Ladders Strongly Coupled by Zener Resonances in Semiconductor Superlat- 
tices, Phys. Rev. Lett. 82 (1999) 3120 
121B. Rosam, D. Meinhold, F. L0ser, V. G. Lyssenko, S. Glutsch, F. Bechstedt, F. Rossi, K. 

Kohler, and K. Leo, Field-Induced Delocalization and Zener Breaakdown in Semiconductor 

Superlattices, Phys. Rev. Lett. 86 (2001) 1307. 
122.B. P. Anderson, T. L. Gustavson, and M. A. Kasevich, Atom Trapping in Nondissipative 

Optical Lattices, Phys. Rev. A 53 (1996) R3727 
123.S. R. Wilkinson, C. F. Bharucha, K. W. Madison, Qian Niu, and M. G. Raizen, Observation 

of Atomic Wannier-Stark Ladders in an Accelerating Optical Potential, Phys. Rev. Lett. 76 

(1996) 4512 

124.C. F. Bharucha, K. W. Madison, P. R. Morrow, S. R. Wikinson, B. Sundaram, and M. G. 
Raizen, Observation of Atomic Tunneling from an Accelerating Optical Potential, Phys. Rev. 



93 



A 55 (1997) R857 

125.K. W. Madison, M. C. Fischer, and M. G. Raizen, Observation of the Wannier- Stark Fan 

and the Fractional Ladder in an Accelerating Optical Lattice, Phys. Rev. A 60 (1999) R1767 
126B. P. Anderson and M. A. Kasevich, Macroscopic Quantum Interference from Atomic Tunnel 

Arrays, Science 282 (1998) 1686 
127.S. Friedel, C. D'Andrea, J. Walz, M. Weitz, and T. W. Hnsch, C02-Laser Optical Lattice 

with Cold Rubidium Atoms, Phys. Rev. A 57 (1998) R20 
128.S. K. Dutta, B. K. Too, and G. Raithel, Tunneling Dynamics and Gauge Potentials in Optical 

Lattices, Phys. Rev. Lett. 83 (1999) 1934 
129.M. B. Dahan, E. Peik, J. Reichel, Y.Castin, and C. Salomon, Bloch Oscillations of Atoms 

in an Optical Potential, Phys. Rev. Lett. 76 (1996) 4508 
130Ij-Guidoni and P. Verkerk, Direct Observation of Atomic Localization in Optical Superlat- 

tices, Phys. Rev. A 57 (1998) R1501 

131. K. Vant, G. Bah, H. Amman, and N. Christcnsen, Experimental Evidence for the Role of 
Cantori as Barriers in Quantum Systems, Phys. Rev. E 59 (1999) 2846 

132. T. Miiller-Scydhtz, M. Hartl, B. Brezger, H. Hnsel, C. Keller, A. Schnetz, R. J. C. Sprccnw, 
T. Pfau, and J. Mlynck, Atoms in the Lowest Motional Band of a Three- Dimensional Optical 
Lattice, Phys. Rev. Lett. 78 (1997) 1038 

133. C. S. Adams and M. Sigel and J. Mlynck, Atom Optics, Phys. Rep. 240 (1994) 143 

134. H. Wallis, Quantum Theory of Atomic Motion in Laser Light, Phys. Rep. 255 (1995) 203 
135.S. Bernet, R. Abfalterer, C. Keller, M. K. Obcrthaler, J. Schmiedmayer, and A. Zeilinger, 

Matter Waves in Time- Modulated Complex Light Potentials, Phys. Rev. A 62 (2000) 023606 

136. Qian Niu, Xian-Geng Zhao, G. A. Georgakis, and M. G. Raizen, Atomic Landau-Zener 
Tunneling and Wannier-Stark ladders in Optical Potentials, Phys. Rev. Lett. 76 (1996) 
4504 

137. M. G. Raizen, C. Salomon, and Qian Niu, New light on quantum transport. Physics Today 
July (1997) 30 

138. K. W. Madison, M. C. Fisher, R. B. Diener, Qian Niu, and M. G. Raizen, Dynamical Bloch 
Band Supression in an Optical Lattice, Phys. Rev. Lett. 81 (1998) 5093 

139. G. Monsivais, M. del Castillo- Mussot, and F. Claro, Stark-Ladder Resonances in the Prop- 
agation of Electromagnetic Waves, Phys. Rev. Lett. 64 (1990) 1433 

140. C. M. de Sterke, J. N. Bright, P. A. Krug, and T. E. Hammon, Observation of an optical 
Stark ladder, Phys. Rev. E 57 (1998) 2365 

141. U. Peschel, T. Pertsch, and F. Lederer, Optical Bloch Oscillations in waveguide arrays. Opt. 
Lett. 23 (1998) 1701 

142. G. Lenz, L Talanina, and C. Martijin de Sterke, Bloch Oscillations in an Array of Curved 
Optical Waveguides, Phys. Rev. Lett. 83 (1999) 963 

143A. Kavokin, G. Malpuech, A. Di Carlo, P. Lugli, and F. Rossi, Photonic Bloch Oscillations 

in laterally confined Bragg Mirrors, Phys. Rev. B 61 (2000) 4413 
14411. Morandotti, U. Peschel, J. S. Aitchinson, H. S. Eisenberg, and Y. Silberberg, Experimental 

Observation of Linear and Nonlinear Optical Bloch Oscillations, Phys. Rev. Lett. 83 (1999) 

4756 

145.T. Pertsch, P. Dannberg, W. Elflein, A. Bruer, and F. Lederer, Optical Bloch Oscillations 
in Temperature Tuned Waveguide Arrays, Phys. Rev. Lett. 83 (1999) 4752 



94 



146. V. I. Sankin, I. A. Stolichnov, and A. A. Mal'tsev, Strong Wannier-Stark localization effects 
in 6H and 4H silicon carbide polytypes, Pis'ma Zh. Tekh. Fiz. 22 (1996) 881 

147. V. I. Sankin, I. A. Stolichnov, Negative differential conduction in the Block oscillations regime 
in the hexagonal silicon carbide polytypes 4-H, 6H, and 8H, Superlattices and Microstructures, 
23 (1997) 999 

148. J. L. Mateos, G. Monsivais, Stark-ladder resonances in elastic waves, Physica A 207 (1994) 
445 

149M. Gliick, A. R. Kolovsky, H. J. Korsch, and N. Moiscyev, Calculation of Wannier-Bloch 
and Wannier-Stark states, Eur. Phys. J. D 4 (1998) 239 

150. M. Gliick, A. R. Kolovsky, and H. J. Korsch, Block particle in presence of dc and ac fields, 
Phys. Lett. A 249 (1998) 483 

151. M. Gliick. A. R. Kolovsky, and H. J. Korsch, Chaotic Wannier-Bloch resonance states, Phys. 
Rev. E 58 (1998) 6835 

152. M. Gliick, A. R. Kolovsky, and H. J. Korsch, A truncated shift- operator technique for the 
calculation of resonances in Stark systems, J. Phys. A 32 (1999) L49 

153M. Gliick, A. R. Kolovsky, and H. J. Korsch, Lifetime statistics for a Bloch particle in ac 

and dc fields, Phys. Rev. E 60 (1999) 247 
154.M. Gliick, A. R. Kolovsky, and H. J. Korsch, Bloch particle in presence of dc and ac fields: 

Statistics of the Wigner delay time, Phys. Rev. Lett. 82 (1999) 1534 
155M. Gliick, A. R. Kolovsky, and H. J. Korsch, Lifetime of Wannier-Stark states, Phys. Rev. 

Lett. 83 (1999) 891 

156.M. Gliick, A. R. Kolovsky, and H. J. Korsch, Perturbation theory for Wannier resonance 

states affected by ac- field, Phys. Lett. A 258 (1999) 383 
157M. Gliick, M. Hankel, A. R. Kolovsky, and H. J. Korsch, Wannier-Stark Ladders in Driven 

Optical Lattices, Phys. Rev. A 61 (2000) 061402(R) 
158.M. Gliick, A. R. Kolovsky, and H. J. Korsch, Fractal stabilization of Wannier-Stark reso- 
nances, Europhys. Lett. 51 (2000) 255 
159M. Gliick, M. Hankel, A. R. Kolovsky, and H. J. Korsch, Induced transitions between Wannier 

ladders, Journal of Optics B: Quantum and Scmiclassical Optics 2 (2000) 612 
160.M. Gliick, A. R. Kolovsky, and H. J. Korsch, Resonant tunneling of Wannier- Stark- states. 

Journal of Optics B: Quantum and Semiclassical Optics 2 (2000) 694 
161 M. Gliick, A. R. Kolovsky, and H. J. Korsch, A quantum cable car for Wannier-Stark ladders, 

Phys. Lett. A 276 (2000) 167 
162.M. Gliick, A. R. Kolovsky, and H. J. Korsch, About Universality of Lifetime Statistics in 

Quantum Chaotic Scattering, Physica E 9 (2001) 478 
163M. Gliick, A. R. Kolovsky, H. J. Korsch, and F. Zimmer, Wannier-Stark resonances in 

semiconductor superlattices, Phys. Rev. B (2002), to appear. 
164.M. Gliick, Wannier-Stark resonances, PhD thesis, Universitat Kaiserslautern, 2000. 
165M. Abramowitz and L A. Stegun, Handbook of Mathematical Fkmctions, Dover Publications, 

Inc., New York, New York, 1972 

166. N. Ashby and S. C. Miller, Electric and Magnetic Translation Group, Phys. Rev. 139 (1965) 
A428 

167. Qian Niu, Effect of an electric field on a split Bloch band, Phys. Rev. B 40 (1989) 3625. 
168A. J. F. Siegert, On the derivation of the Dispersion Formula for Nuclear Reactions, Phys. 



95 



Rev. 56 (1939) 750 

169. M. Wagner and H. Mizuta, Complex- energy analysis of intrinsic lifetimes of resonances in 
biased multiple quantum wells, Phys. Rev. B 48 (1993) 14393 

170. E. Hernandez, A Jauregui, and A. Mondragon, Degeneracy of Resonances in a Double Bar- 
rier Potential, J. Phys. A 33 (2000) 4507 

171. M. Philipp, P. von Brentano, G. Pascovici, and A. Richter, Frequency and Width Crossing 
of Two Interacting Resonances in a Microwave Cavity, Phys. Rev. E 62 (2000) 1922 

172. U. Fano, Effects of configuration interaction on intensities and phase shifts, Phys. Rev. 124 
(1961) 1866 

173. H. A. Kramers,, in R. Stoops, editor, Les Particules Elementaires, Proceedings of the Eighth 
Solvay Conference, New York, 1950. Wiley 

174. H. A. Kramers, Collected Scientific Papers, North-Hohand, Amsterdam, 1956 

175. W. C. Hcnncbcrger, Perturbation Method for Atoms in Intense Light Beams, Phys. Rev. 
Lett. 21 (1968) 838 

176. C. K. Choi, W. C. Henneberger, and F. C. Sanders, Intensity- dependent Ionization Potentials 
for H and He in Intense Laser Beams, Phys. Rev. A 9 (1974) 1895 

177. J. Zak, Finite Translations in Time and Energy, Phys. Rev. Lett. 71 (1993) 2623 

178J. Zak, Quasienergy states for a Bloch electron in a constant electric field, J. Phys. Cond. 
Mat. 8 (1996) 8295 

179.F. G. Bass and A. P. Tetervov, Hight-frequency phenomena in semiconductor superlattices, 

Phys. Rep. 140 (1986) 237 
180D. R. Hofstadter, Energy Levels and Wave Functions of Bloch Electrons in Rational and 

Irrational Magnetic fields, Phys. Rev. B 14 (1076) 2239 

181. F. Bensch, H.-J. Korsch, and N. Moiseyev, Simple method for constructing the ionization 
spectra of driven time-periodic Hamiltonians, Phys. Rev. A 43 (1991) 5145 

182. F. T. Smith, Lifetime Matrix in Collision Theory, Phys. Rev. B 118 (1960) 349 

183. W van Dijk, F. Kataoka and Y. Nogami, Space-time evolution of a decaying quantum state, 
J. Phys. A 32 (1999) 6347 

184. W van Dijk and Y. Nogami, Novel Expression for the Wave Function of a Decaying Quantum 
System, Phys. Rev. Lett. 83 (1999) 2867 

185. J. R. Taylor, Scattering theory, John Wiley, New York, 1972 

186. A. S. Parkins and D. F. Wahs, The Physics of Trapped Dilute-Gas Bose-Einstein Conden- 
sates, Phys. Rep. 303 (1998) 1 

187. M. L. Chiofalo and M. P. Tosi, Output from Bose Condensates in Tunnel Arrays: The Role 
of Mean-Field Interactions and of Transverse confinement, Phys. Lett. A 268 (2000) 406 

188. K. Berg-S0rensen and K. M0lmer, Bose-Einstein Condensates in Spatially Periodic Poten- 
tials, Phys. Rev. A 58 (1998) 1480 

189. Dae-Il Choi and Qian Niu, Bose-Einstein Condensates in an Optical Lattice, Phys. Rev. 
Lett. 82 (1999) 2022 

190. Biao Wu and Qian Niu, Nonlinear Landau-Zener Tunneling, Phys. Rev. A 61 (2000) 023402 

191. J. Javanainen, Phonon Approach to an Array of Traps Containing Bose-Einstein Conden- 
sates, Phys. Rev. A 60 (1999) 4902 

192.0. Zobay and B. M. Garraway, Time- dependent Tunneling of Bose-Einstein Condensates, 
Phys. Rev. A 61 (2000) 033603 



96 



193. M. M. Ccrimcle, M. L. Chiofalo, F. Pistella, S. Succi, and M. P. Tosi, Numerical Solution of 
the Gross-Pitaevskii Equation Using an Explicit Finite-Difference Scheme: An Application 
to Trapped Bose-Einstein Condensates, Phys. Rev. E 62 (2000) 1382 

194. F. Haake, Quantum Signatures of Chaos, Springer, New York, 1991 

195. T. Guhr, A. Miiller-Groeling, and H. A. Weidenmiiller, Random-Matrix Theories in Quan- 
tum Physics: Common Concepts, Phys. Rep. 299 (1998) 189 

196. C. W. J. Beenakker, Random-Matrix Theory of Quantum Transport, Rev. Mod. Phys. 69 

(1997) 731 

197.0. Bohigas, M. J. Giannoni, and C. Schmit, Characterization of Chaotic Spectra and Uni- 
versality of Level Fluctuation Laws, Phys. Rev. Lett. 52 (1984) 1 

198. N. Lehmann, D. V. Savin, V. V. Sokolov, and H.-J. Sommers, Time Delay Correlations in 
Chaotic Scattering: Random Matrix Approach, Physica D 86 (1995) 572 

199. Y. V. Fyodorov and H.-J. Sommers, Statistics of resonance poles, phase shifts and tim,e 
delay in quantum chaotic scattering: Random matrix approach for systems with broken time- 
reversal invariance, J. Math. Phys. 38 (1997) 1918 

200. P. Scba, Random Matrix Theory and Mesoscopic Fluctuations, Phys. Rev. B 53 (1996) 13024 

201. P. Scba, K. Zyczkowski, and J. Zakrewski, Statistical Properties of Random Scattering Ma- 
trices, Phys. Rev. E 54 (1996) 2438 

202. H.-J. Sommers, Y. V. Fyodorov, and M. Titov, S-Matrix Poles for Chaotic Quantum Systems 
as Eigenvalues of Complex Symmetric Random Matrices: From Isolated to Overlapping 
Resonances, J. Phys. A 32 (1999) L77 

203. H. Ishio and J. Burgdorfer, Quantum Conductance Fluctuations and Classical Short-Path 
Dynamics, Phys. Rev. B 51 (1995) 2013 

204. L. Wirtz, J.-Z. Tang, and J. Burgdorfer, Geometry-Dependent Scattering through Ballistic 
Micro structures: Semiclassical Theory beyond the Stationary-Phase Approximation, Phys. 
Rev. B 56 (1997) 7589 

205. H. Ishio, Resonance Poles and Width Distribution for Time-Reversal Transport through 
Mesoscopic Open Billards, Phys. Rev. E 62 (2000) R3035 

206. R. Bliimel, Existence of a Ericson Regime in Stretched Helium, Phys. Rev. A 54 (1996) 5420 

207. V. A. Mandelshtam and H. S. Taylor, The quantum resonance spectrum of the molecular 
ion for J = 0. An accurate calculation using fUter-diagonalization, J. Chem. Soc, Faraday 
Trans. 93 (1997) 847 

208. V. A. Mandelshtam and H. S. Taylor, Spectral Analysis of Time Correlation Function for a 
Dissipative Dynamical System Using Filter Diagonalization: Application to Calculation of 
Unimolecular Decay Rates, Phys. Rev. Lett. 78 (1997) 3274 

209P. Borgonovi, L Guarneri, and D. L. Shepelyansky, Statistics of Quantum Lifetimes in a 
Classically Chaotic System, Phys. Rev. A 43 (1991) 4517 

210. G. Casati, G. Maspero, and D. L. Shepelyansky, Relaxation Process in a Regime of Quantum 
Chaos, Phys. Rev. E 56 (1997) R6233 

211. G. Casati, L Guarneri, and G. Maspero, Fractal Survival Probability Fluctuations, Phys. 
Rev. Lett. 84 (2000) 63 

212. T. Kottos and U. Smilansky, Quantum Chaos on Graphs, Phys. Rev. Lett. 79 (1997) 4794 

213. T. Kottos and U. Smilansky, Chaotic Scattering on Graphs, Phys. Rev. Lett. 85 (2000) 968 
214A. J. Lichtenberg, M. A. Lieberman, Regular and Chaotic Dynamics, (Springer, Berlin, 1983) 



97 



215.R. Graham, M. Schlautmann, and P. Zoller, Dynamical localization of atomic-beam deflection 

by a modulated standing light wave, Phys. Rev. A 45 (1992) R19 
216JVI. Pozniak, K. Zyczkowski, and M. Kus, Composed ensembles of random unitary matrices, 

J. Phys. A 31 (1998) 1059 
217.K. Zyczkowski and H.-J. Sommers, Truncation of Random Unitary Matrices, J. Phys. A 33 

(2000) 2045 

218D. V. Savin, Y. V. Fyodorov, and H.-J. Sommers, Reducing Nonideal to Ideal Coupling in 
Random Matrix Description of Chaotic Scattering: Application to the Time-Delay Problem, 
Phys. Rev. E 63 (2001) 035202(R) 

219. D. V. Savin, private communication 

220. D. V. Savin and V. V. Sokolov, Quantum versus Classical Decay Laws in Open Chaotic 
Systems, Phys. Rev. E 56 (1997) R4911 

221. J. C. Robinson, C. Bharucha, F. L. Moore, R. Jahnkc, G. A. Gcorgakis, Q. Niu, and M. 
G. Raizen, Study of Quantum Dynamics in the Transition from Classical Stability to Chaos, 
Phys. Rev. Lett. 74 (1995) 3963 

222. A. R. Kolovsky, Bloch oscillations of the atoms in a standing laser wave, unpubhshed 

223. M. Gliick, F. Keck, A. R. Kolovsky, and H. J. Korsch, Wannier-Stark states of a quantum 
particle in 2D lattices, Phys. Rev. Lett. 86 (2001) 3116 

224.S. Flach, O. Yevtushenko, and Y. Zolotaryuk, Directed current due to broken time-space 

symmetry, Phys. Rev. Lett. 84 (2000) 2358 
225.T. Dittrich, R. Ketzmerick, M. F. Otto, and H. Schanz, Classical and quantum transport in 

deterministic Hamiltonian ratchets, Ann. Phys. (Leipzig) 9 (2000) 755 (2000) 
226H. Schanz, M.-F. Otto, R. Ketzmerick, and T. Dittrich, Classical and quantum Hamiltonian 

ratchets Phys. Rev. Lett. 87 (2001) 070601 
227.1. Goychuk and P. Hanggi, Directed Current without Dissipation: Reincarnation of a 

Maxwell-Loschmidt Demon, Lecture Notes in Physics: "Stochastic Processes in Physics", 

J. Freund and T. Poschl, Eds., Vol. 557, pp. 7-20, Springer Verlag, Berlin, Heidelberg, 2000 
2281. Goychuk and P. Hanggi, Minimal Quantum Brownian Rectifiers, J. Phys. Chem.105 (2001) 

6642 

229JI. Utermann, T. Dittrich, and P. Hanggi, Tunneling and the Onset of Chaos in a Driven 

Bistable System, Phys. Rev. E 49 (1994) 273 
230.S. Tomsovic and D. UUmo, Chaos assisted tunneling, Phys. Rev. E 50 (1994) 145 

231. V. Averbukh, N. Moiseyev, B. Mirbach and H. J. Korsch, Dynamical Tunneling through a 
Chaotic Region - A Continuously Driven Rigid Rotor, Z. Phys. D 35 (1995) 247 

232. W. K. Hensinger, H. Haffner, A. Browaeys, N. R. Heckenberg, K. Helmerson, C. McKen- 
zie, G. J. Milburn, W. D. Phillips, S. L. Rolston, H. Rubinsztein-Dunlop, and B. Upcroft, 
Dynamical tunnelling of ultracold atom,s, Nature 412 (2001) 52 

233D. A. Steck, W. H. Oskay, and M. G. Raizen, Observation of Chaos-Assisted Tunneling 

between Islands of Stability, Science 293 (2001) 274 
234B. Goss Levi, Atoms hop between islands of regular motion in a sea of chaos. Physics Today, 

August 2001, 15 



